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ABSTRACT 


An investigation of neutron neutron scattering 
has been made in the framework of a kinematically complete 
D(n,2n)p experiment at 21.5 MeV incident energy where the 
three body breakup reaction was observed under n-n Quasi 
Free Scattering geometry. 

To improve the low counting rate problem involved 
in the neutron induced reaction a cryogenic facility, 


/ 23 MeV neutrons/sr sec, was 


capable of producing 4 x 10 
designed and comissioned in this project. A more effective 
use of the neutron flux produced was achieved by exposing 
to it an array of deuterium targets in the form of NE230 
deuterated scintillators. A set of three such scintillators 
provided simultaneously, data at three different scattering 
geometries. Data at a total of nine different quasi-free 
scattering angle pairs have been obtained in three 
experimental runs. 

The proton pulse in each target scintillator was 
used as a reference timing signal in the double time-of- 
flight spectrometer adding to a total of eight parameters 
processed in coincidence by the electronics for each event. 
Events from the different NE230 scintillators were sorted 


by means of four-way routing and recorded event-by-event 


On a magnetic tape. 
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The experimental results are reproduced reasonably 
well by a separable potential approximation of the exact 
three-body theory in contrast with up to 80% discrepancies 


reported earlier in the p-p quasi-free scattering experi- 


ments at this energy. 
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The left and right monitor spectra 
corresponding to the entire runs #1, 
#2, and #3. The inset is showing 
the contribution of n-C elastic 
scattering into the n-D elastic 
monitor peak and the reconstruction 
of their superposition. 


The neutron detector efficiency curves 


used ine this work. 


Comparison of experimental and 


theoretical differential cross sections 
for the reaction n,td > n,+no+p in the 


projection onto Eh axis wn. Rune 1) 
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The dashed line indicates the proton 

energy (right scale). 


Comparison of experimental and 


theoretical differential cross sections 


(solid line) for the reaction 
njtd ae n,tnotp in projection onto the 


kinematical locus length S for the 
three angle pairs of Run #2. 


Comparison of experimental and 


theoretical differential cross sections 


(solid line) for the reaction 
njtd a n,tnotp in projection onto the 


kinematical locus length S for the 
three angle pairs of Run #3. 
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Comparison of relevant energies and 
both experimental and theoretical 
differential cross sections for the 
reaction njtd all yato  P at Ejne2l.9 MeV 
with 4 Ors >5=180 8,=8,=25 F 
(Upper) Absolute experimental data 
(step line) and theoretical curves. 
(Lower) Relative energies a p? 
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As melHEORY 


Avvler THEORY OF NUCLEAR: FORGES 
A.1.1 Meson Exchange Theory of the Nucleon Nucleon Force 
The modern theory of nuclear forces explains the 
nucleon-nucleon nuclear or "strong" interaction in terms of 
exchange of mesons. A simplified illustration of this 
hypothesis can be given on the basis of Heisenberg's rela- 
tions of uncertainty in the Quantum Mechanical micro- 
dimensions (He-27). Accordingly, the total energy E of an 
isolated nucleon is allowed to fluctuate with an uncertainty 


AE related to its time duration Tt by 
AEs t = ti (A-1) 
A Pas Pe lanck. Ss Cons bane cm 


AE can be expressed in terms of a rest energy mee 
of a particle, meson, which the nucleon can constantly emit 
and reabsorb within time intervals t. The meson is 


"allowed" to travel a distance 
- e ~ —E Vv -_ 
Reape hy ie : (A-2) 


during such interval if v is its speed. Thus the nucleon 
could be viewed as a "gallaxy" of mesons, the heavy hadrons 


concentrated in its centre core and the lightest mesons 
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thinning off at the longer ranges R from the center. If two 
nucleons move to within the meson emission range of each 
other a meson emitted by one nucleon can be absorbed in the 
other. The corresponding change in the state of each indivi- 
dual nucleon due to the meson exchange must then be described 
as their interaction in a two-nucleon system to conserve 
energy. 

Li Vic in equation (A-2) . for orientation, one 


obtains: 


= i | (A-3) 


The range of nuclear interaction is thus equal to the Compton 
wave-length Xe of the exchange particle. For R = taut 1.4F 
(the observed range of nuclear forces) the corresponding rest 


mass of the exchanged meson is m = 270 m » which is 


electron 
the known II meson. 

A number of mesons of different masses which can 
be exchanged in the nucleon-nucleon interaction are known 
today (7, n, 9p, w, >). (see Table 1) 

Depending on the proximity of the two interacting 
nucleons one pion, two pions, n, 9p, w, » mesons and their 
combinations can be exchanged in the interaction making its 


mathematical expression (Pa-70, Lo-68) correspondingly 


complex if not impossible. 
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particle spin 

orbital angular momentum 
particle isospin 

parity G = 
ye 


particle parity P = -(-1 


charge conjugation parity eigenvalue 


(Re-76) for more information 
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The simplest of the exchange forces, the One Pion 
Exchange Potential (OPEP), is based on the derivation of 
Yukawa (Yu-35) who first hypothesized the meson exchange 


theory of nuclear forces. Yukawa's potential 
make 


U(r) = -g° & (A~4) 


(where g is the strong pion-nucleon coupling constant and 
] 


zs “1s. the Compton wavelength of pion), assuming the exchange 
of one neutral pion, is applicablhe omly for so called long 
range interaction, i.e., an interaction in which the two 
nucleons do not come closer than = 1.5F. At closer ranges 
combinations of simultaneous meson exchange modes alter the 
interaction drastically. The pseudoscalar or vector 
Character of mesons in spin-isospin space included in the 
more recent theories is able to account for some observed 
features of nuclear forces, which Yukawa's central potential 
fails to explain, (e.g., the spin-dependence of nuclear 
forces their non-central character, the repulsive core of 
the potential, etc.). One of the most sensitive probes into 


the details of the meson exchange theory is to test the 


charge independence and charge symmetry of nuclear forces. 


A.1.2 Principle of Charge Independence and Charge Symmetry 
It is known that the hadronic part of the nucleon- 
nucleon force is not only charge symmetric but also charge 


independent (He-69). It is also known that these postulated 
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symmetries are proper to the strong or hadronic interaction 
alone. They are both broken by the electromagnetic and weak 
interactions. The strength of electromagnetic interactions 
is characterized by the fine structure constant o = 1/137. 
As, on the grounds of the meson exchange theory, the hadronic 
part in the nucleon-nucleon interaction is always accompanied 
by long range or short range electromagnetic interactions, 
it is expected that, in general, the charge independence and 
charge symmetry of the nucleon-nucleon forces will be broken 
by a few percent. 

Charge independence applied to the nucleon implies 
that once Coulomb and other small electromagnetic effects 
are removed, the n-n, n-p and p-p forces are equal in the 
Sameuc la vest 

Charge symmetry is weaker implying equality only 
of n-n and p-p forces. 

The breaking of charge independence is primarily 
due to the electromagnetic mass differences of the exchanged 


© ~see TABLE 1). These effects do not 


mesons (e.g., ae T 
give rise to charge symmetry breaking. 

Although there is a clear evidence for the 
violation of charge independence (He-69), this is not so 
for charge symmetry. On theoretical grounds, charge- 
Symmetry violating effects must occur for nuclear forces, 


unless there are accidental cancellations. Such symmetry- 


breaking effects follow from the theoretically predicted 
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mixing of isoscalar and isovector mesons (Go-/70). 

There is indirect evidence for charge symmetry 
breaking (O0k-67, Ne-71) in low-energy nuclear physics. This 
comes primarily from the measured binding-energy differences 


of the mirror nuclei 3e- 34 and Tee G 


a. The theoretical 
analyses require an extra attraction of the low energy n-n 
system relative to the p-p system of roughly 1/2%. 

Attempts to test the validity of charge symmetry 


through direct study of the low energy nucleon-nucleon 


System are described in the next chapter. 


A.1.3 The Low Energy NN System 

The comparison of NN bound states and scattering 
cross sections, provides a direct test of charge symmetry 
and charge independence. Unfortunately, only the np system 
has a bound state, the deuteron, which occurs for an 1=0 
spin triplet. At low energies (lab. kinetic energy < 10 
MeV) only s-state (2=0) NN scattering is important and 
contributions from higher partial waves can be neglected. 
The corresponding spin-isospin states allowed by the Pauli 
principle are the spin singlet ‘Sj7 isospin triplet (T=1) 
in the nn, np and pp scattering plus the spin triplet oe 
isospin singlet (T=0) scattering in np system. The latter 
contribution can be easily removed from the np scattering 
amplitude with the help of slow-neutron coherent scattering 


in Hydrogen (Wi-63). Thus at low energies the cross 
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sections for the three T=1 members of NN scattering can be 
compared in the same space-spin states Sa One advantage 
of the low energy scattering is, that the scattering 


cross section can be expressed in terms of only 2=0 phase 


shift on as 
qo 2 sin? 60 (k) (A-5) 
k ... C.M. momentum 


Further advantage is, that this phase shift can be expressed 


at low energies by the expansion 


ig Ko 


1 3,4 
aie ( ah eg) 


: ei 
Ke = ke cot om 5 7 alee: + 


; (A-6) 


The comparison of cross sections is thus reduced at low 
energies to a comparison of just two-parameters, the 
scattering length (a) and the effective range bre 

It is known experimentally that the nn, pn and pp force in 
the spin singlet 4S is attractive and the states are 
almost bound making the scattering length large and negative. 
Under such conditions a slight change AV in the depth of 
nucleon-nucleon potential V reflects strongly magnified 

in the corresponding change Aa of the scattering length as 
the two variables are related by 
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Since the scattering length a is about an order of 
magnitude larger than the effective range ine of the NN 
potential in case of "ge state a large magnification of the 
charge dependent effects is provided by the comparison of 
nn, pp and np scattering lengths. 

Unfortunately, two major hardships beset the 
direct experimental comparison of the three NN scattering 
States and these are responsible for existing uncertainty 
about the EHARE Sy mnethy of NN forces to date. These are: 
(i) the nonexistence of a free-neutron target for direct 

n-n scattering, and 
(ii) the presence of the Coulomb force in p-p scattering. 

There is no long range electromagnetic force 
present in the n-p interaction. Careful analysis of the 
experimental data below 5 MeV, including epithermal neutron 
scattering, binding energy of the deuteron and coherent 


scattering of cold neutrons (No-68) yield 


a 226 7146 4 0.0127 F 
mp } (A-8) 
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rnp 
The p-p scattering experiment, being the easiest 
to perform, has been carried out to a great accuracy. The 


result (He-67) 
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It 


- 7.804 + 0.006 F 
ma } (A-9) 
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was used as a basis for the theoretical prediction of the 


charge-dependent nn scattering length a which incorporated 


nn 
the effects of p-w and n-1 mixing in the meson exchange 
model (He-72). First, however, the presence of the long 
range electromagnetic interaction in result (A-9) has to be 
“turned off", which can only be done theoretically. The 
result of the theoretical subtraction of Coulomb and vacuum 


polarization, which is not model-independent, in the pp 


scattering length by Henley et al. (He-72) is: 


Nese Ase] 71 A-10 
a “fi 7.10 ( ) 


The effective range is insensitive to Henley's correction. 
The uncertainty in result (A-10) was quoted by Henley as 
+ 0.20 F. However, in his recent study, Sauer (Sa-76) 
discloses that, the model dependence of the subtraction of 
the electromagnetic effects is, in fact, extremely serious. 
He concludes that the value of -17 F for the nuclear pp 
scattering length is still most probable, but its model 
dependence arising from different parametrization of the 
wave function at small relative distances is disappointingly 
large. 

On the basis of result (A-10) Henley (He-72) made 
a theoretical prediction of two possible values for the 
charge dependent (i.e., such as should be found in experi- 
ment) neutron-neutron scattering length, depending on the 


sign of the ratio Sick of w and p meson strong coupling 
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constants, which could not be predicted by SU(3) symmetry, 


SN eo HE ma ae Atle hig mac) 
Se (A-11) 
= WG 0p ae g/9 <Q 


He then uses the evidence of an extra attraction needed in 


S 3 


the n-n system over p-p system, to explain the ~He-~H and 


elspa a 


Ca binding energy differences, as an argument for his 
Final choices ot-tne — 1/795 value over — 1520 FE “for: the 


predicted a This value has to be verified experimentally 


aS 
to prove the validity of the meson exchange model used for 
the prediction. Henley's value of the nn effective range 


from the summary of present knowledge (He-69) is 
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Due to nonexistence of a free-neutron target the 
nn scattering parameters can only be obtained indirectly in 
many body (mainly few-particle breakup) reactions where the 
nn subsystem appears in the final state. The results of 
such analyses are dependent on, more or less justified, 
approximations accounting for the effect of the many body 
framework onto the semi-isolated n-n system. The simplest 
of the many-body reactions providing the same framework for 
the comparison of nn, np and pp interactions in its final 
states is the three body nucleon induced breakup of the 


deuteron. 
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A.2 THREE-NUCLEON BREAK-UP OF THE DEUTERON 

In recent years, interest in experimental and 
theoretical investigations of three-nucleon systems, expecially 
the break-up reactions “H(p.2p)n and 2u(n,2n)p, has been 
continuously increasing. The aim of all the studies on this 
subject has been three-fold (S1-70): 

(i) A systematic analysis of the most important reaction 
mechanisms; that is, final-state interactions in two- 
nucleon subsystems (f.s.i.) and quasi-free nucleon- 
nucteonascattering (q.f.5.). 

(ii) A quantitative description of the data by approximate 
solutions of the exact three-body theory. 

(iii)The determination of nucleon-nucleon scattering 
parameters. 

The advantage of the nucleon and deuteron three- 
body reaction is its ability to provide the comparison of 

the three charge states of the nucleon-nucleon Sig scattering 

under equal conditions in the form of nt+td > ntnt+p or 

ptd + ptp+n break-up. The presence of the three particle 

framework in the study of the nucleon-nucleon subsystem is 

not such a disadvantage since the exact three-body theory 

(Fa-61) and its solutions (Lo-64, Eb-72) have been developed. 

On the contrary, the study of the three-body system can reveal 

new information regarding details of the two-particle 


interactions, which cannot be readily observed in the study 
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of the two-body system. In particular the short range 
behaviour of the nuclear force, the off-shell effects and, 
in addition, the possible existence of the three-body 


force (Br-74) can be studied in the three-body system. 


A.2.1 Kinematics of Three-Particle Reactions and the 
Differential Cross Section 


We consider reactions of the following form: 
OR te [9s aes 


The nine scalar momentum components of the final state 
completely specify the kinematics of the system. Four 
components can be eliminated by applying conservation 
of energy and momentum. Thus the measurement of five 
independent scalar variables sufficies. 

In coincidence experiments, two of the particles 
(by our convention, particles 1 and 2) are detected at fixed 
angles. Thus six scalar variables are measured: the beam 
energy, two final-state particle energies, two angles 
with respect to the beam axis, and the azimuthal angle between 
the two detected particles. This constitutes a kinematic 
over specification of the system, which is useful in reducing 
background. Conservation of energy and momentum restricts 
the energies of the two detected particles to values given 


by the following equation: 
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Q ri (] 7 m,/m,)E, + (1 ai m/m,)E, = 
(1 - mo /ms)E. - 2c0os 8, [(m,m /ms ye £, 11/4 


85 (mpm, /m5)E E,1!/* 


i 
PO 
(2) 
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+ 2cos 6, £(m,m,/ms)E,E,1'/? Cnet) 


where cos 815 = cos 6, cos 0, + sin 8, sin 8, cos (Wy -V5) 


The O55 Was m. 


F and E. are the polar angles, azimuthal angles, 


masses, and energies of particle i; Q is the Q-value of the 
reaction. This is the equation of a fourth-degree curve 

in the E,-E. experimental space (indicated as dotted and 
solid curves in the two-dimensional plots of Chapter C). 
Particle 3 is free to assume various momenta consistent with 
each E,> E, SOuUtTON Ot bq.) (A=-13}. "bach point om the 


E, = a) curve corresponds to a definite relative energy 


2 
between each pair of particles. The occurrence of a reso- 
nance at a particular relative energy results in an enhance- 
ment over the normal phase space or a peak in the three-body 
cross section at the appropriate point on the E. ~ F(E,) 
curve. An experiment at fixed angles and fixed incident 
energy may contain several peaks corresponding to resonances 
Gravar Vous Internal relative energies (bids. Oona C5 a5))4 
The intensity distribution I(s ,8) ,85) along the 


curve contains the full physical information. The experimen- 


» ys ~~ 4 ~ nal ~ 7 _ v = @ 
oF ts | ks et: On | - 

. iy % -e ; Le a 
= : . - i > 


; u<. ~ 
As salem! my" 2ieaunt ea 


tates hielo pT 608. 4a(gmiga'= tT)? 


s\i eee ts. 

[A gatemongo lt ce 2008 
(EF-A) S\ Ifa at Savqmipin)d 7 2098 + 
(g¥ ;W) 209 8 nite (8 ne # <6 203 di 209 = 9,8 209 sien 


-esipns isdtumtss “2afous aplog ond S16 ‘os bie ue i +8 oaT = 
eit *0 aulev-0 ant 2h 0 ft sloiiteg to 2stpvens bas -esceom oe 
aviv saipasb-diyvoh 6 to nOfisuns snt ef cid?) =. noheoesy 

bis bettob 26 bsisarbat) s28ge fednomtreqxs c3- 3 eudoat 

.(2 vetqsdd to esolg fsnorenomtb-ows sax ar zevawa biloe r 

d¢iw toeteatenos sinemom evotrev smbe2s oF sett ef e@ sfotsisd 
adi qo tntoq dos3 .(ET-A) .p3 a nottuloe <3 ed foes 

ypisns aviteter 9tiniteb 6 o2 2bnoq2a¥109 avi ( rit 7) 
-029% 6 To S9Nn9TIUDI0 SAT .asfotiisg +o 1isq tose. ds 
“sonetins no nf esivesy ypyons ov tie l9n wafuotiysg 6 
ybod-s9ynt oft nt 459q & 10 SD6Ge S25 ismvon add ala 
(pat "A 2 adit no tatog stetyqorqqs afd el i 22 
dwabtant bexi? bos 2ofpns boxtt ts peer nA : 
250n6n0Z97 D eNTh ange SAAD: cAnog, Fenste ovate: nite 
(2.2.9 8-9 2077) eaiotane pein , 


14 


tal data are optimally analysed in terms of the arc length 
s(E,,E,) along this curve. The differential cross section is 
expressed by 

5 


ds 2, Q 


Shima (A-14) 
0 


: fi 


with the projectile velocity vies the transition amplitude 


Te, and the phase space factor 
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A.2.2 Transition Amplitude Te, - The Exact Three Body 
Formalism 


Exact equations for the three-body transition 
amplitude were first given by Faddeev (Fa-61). He expressed 
the three-body transition operator formally as a multiple 
scattering series 
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ee 


3 terms 6 terms 12 terms 


Here t is a two-particle transition operator and the two- 
particle transition amplitude from a state of two particles 
with a relative momentum p to a state with relative momentum 
p' would be <p'|t(e)|p>. Here e€ is the two-particle relative 
energy. If the two-particle bound state or a resonance exist 
the transition amplitude can be divided into a separable 

part (2°05 which contains the pole at the bound state 

energy € = - ee or a sharp peak at e€ = 0, and the rest 


without any pole behaviour as: 
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Gpr jareligtipje + tres't{ praipyel) (A-18) 


Here g{p) are form factors and t(c) is a function carrying 
the pole behaviour in the relative energy of the two particles 


€. This can be diagrammatically symbolized as: 


= «a Di c= = Olea pal. (A-18a) 
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For a separable potential of Yamaguchi (Ya-54) 
type v(p',p) = Ag(p')g(p), the remainder of Eq. (A-18) is 
identically zero. The Yamaguchi form factors are connected 


to the scattering length a and the effective range e by 


oO 
-— 
AE} 
~ 
It 
me 
= 
— 
qe) 
ae 
WD 
_— 


(A-19) 


Solutions of the Faddeev equations (A-17) with local 
nucleon-nucleon potentials are not available for the break- 
up reaction, whereas a separable potential expansion (Lo-64), 
especially in the reformulation of the Faddeev equations by 
Sandhas et al., (Sa-67) has led to numerical calculations of 
thepabSouluce dittherential cross. Section, |Ca-7 1, as Bbs72). 

The soluable coupled integral equations for the 
breakup transition amplitude in the Sandhas' rearrangement, 
for separable potentials only, can be expressed by the 


following diagram: 
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SUNN eae seeser OB and canara sie are defined by (A-18a) 
and symbolizes the breakup transition 
amplitude. 
The solution of integral "equation" (A-20) by 

Ebenhoh (Eb-72), the theoretical predictions of which were 
used in this work, has been obtained with the following 
important restrictions to the three-particle interaction: 
(i) Only two-particle S-wave interactions in triplet 

and singlet states of the two-nucleon subsystems 

are. taken into account. 
(ii) The Coulomb repulsion of the p-p subsystem is 

not included (it, however, does not appear in the 

d(n,2n)p reaction). 
(iii) Separable potentials of Yamaguchi type fitted to 

low-energy effective range parameters are used to 


solve the coupled integral equations numerically. 


A.2.3 Previous Studies of Deuteron Break Up Reactions 

So far, experimental studies of three-nucleon 
breakup reactions have shown enhanced differential cross 
sections for certain kinematical situations only, namely 
for a small relative energy ENN of any two nucleons (Final 
State Interaction) and for a minimum momentum transfer to 
one nucleon bound in the deuteron (Quasi-Free Scattering). 


Most of the experimental studies of deuteron breakup reac- 
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tions performed up to date concentrated on these semi two- 
body effects under kinematical conditions, where they would 


dominate the breakup and give enhanced cross sections. 


Meeacel Final State interactions 

Under kinematical conditions where two outgoing 
nucleons from the breakup have small relative energy while 
their energy with respect to the third particle is large, or 
if the third particle cannot effect the interaction of the 
first two particles in the final state, the breakup transition 
amplitude is fully dominated by the two-particle interaction. 
In the approximation of the Watson-Migdal model (Wa-52, Mi-55) 
the probability of all but one two-particle final state is 
neglected and the transition amplitude is expressed by the 


sequence (e.g.} in n + d +n +n + p) 


o=5 aa 


where the production amplitude i is assumed constant for 


fixed angles of detection and ue has a pole behaviour at the 


n 


Pant relative momentum. This reflects in factorization of 


the cross section as 
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is an enhancement factor. As a f.s.i. dominates the breakup 
process very strongly under appropriate kinematic conditions 
there is relatively little deviation between the values of 
the scattering length a, as deduced from experiments by 
"exact" calculations and those obtained with the W.M. 
approximation (Bo-69). While the experiments reproduce the 
established value of anp and the Coulomb corrected Anny? 
there existed large discrepancies in the earlier, so called 
kinematically, incomplete experiments in the value obtained 
ror aa (Sh-73 references in there). These experiments did 
not collect enough kinematical variables for a complete 
definition of the breakup. Instead, they detected one 
particle, analysing a forward direction peak in its angular 
distribution. This corresponded to the case where all 
available kinetic energy is transferred to this particle, 
leaving the remaining two with zero absolute and thus also 
zero relative energy. Approximate solution methods, like 
Born Approximation and W.M. model, used in the analysis of 
these experiments before the advent of the "exact" calcula- 
tions, proved inadequate to cope with the interference of many 
effects involved in such reactions. The a values obtained 
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The results of more recent, kinematically complete, 
experiments (Eg., Ze-70, 72) cluster in the range between 
CECE 14.5 F to - 18.0 F. The summary of these experiments 
| by Vranic et al. (Vr-74) assigns the nn effective range 


parameter values to a Shia wl Oot ateD 29 at oe 


xa Sveqlee wn soeige 


nn 
The kinematically complete experiment which is best 
satisfying the assumptions of the W.-M. approximation is the 
“eaGtiOnan stedeet- cyet—2n.,sThisvhaspbeen. studied by 
Haddock et al. (Ha-65). The result, corrected by Nygren 
(Ny-68), was eer as Len42st41253 Feo (heganalyses,oft 

f.s.i. experiments appear to be rather insensitive to the 


effective range r and its value is usually taken as equal 


nn 


to "op (oak iiecaloulatioans . 


A.2.3.2 Quasi free scattering in the deuteron break up 

At high incident nucleon kinetic energies (> 100 
MeV), where the Compton wavelength of the bombarding nucleon 
is small enough to resolve the two nucleons bound in the 
deuteron, quasi-free scattering can simply be explained as 
almost-free scattering of two of the three nucleons involved 
in the deuteron break up, under kinematic conditions where 
the "spectator" particle is left without any share of the 
available kinetic energy. This explanation, however, is not 
sO applicable at low energies (~ 10-30 MeV) where the deuteron 
binding energy of 2.225 MeV is not negligible. The effect 


spreads into a coherent contribution of a number of coupled 
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multiple scattering terms in the three-body transition matrix 
which involve interactions of all the three nucleons. In 

any case, a broad enhancement peak is observed in the low 
energy deuteron breakup at the geometry, which corresponds 

to the quasi-free scattering, and the process retains 
predominance of the on-shell two-body scattering. 

The semi-two-nucleon kinematics, perturbed by the 
binding of the spectator nucleon in the deuteron, changes 
the laboratory angle between the scattered nucleons from 
~ 90° (which would correspond to a free two-nucleon 
scattering) to a slightly smaller angle, depending on the 
ratio of the C.M. energy to the deuteron binding energy. 
Typically, the geometry of the quasi-free scattering at 
~ 20 MeV incident nucleon energy in the laboratory system 
is 0, Seto % 8. = - 41.5° for the two outgoing nucleon 
scattering angles in the laboratory system. 

The multiple scattering contributions involved 
in the low energy q.f.s. extend the dominance of this effect 
both onto a large range of scattering angles around the pure 
q.f.s. geometry as well as onto a large region, along the 
kinematical locus E2=f(E1), corresponding to one particular 
d.tys<s geometry. 

The study of 9.f.s. has mostly been done for the 
D(p,2p)n reaction in a broad range of bombarding energies 


(4752 '60 MeV; "Va-71; Du=771; 'Pe-70). ‘The “broad peak ‘due “to 
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the p-p q.f.s. observed in these experiments is reproduced, 
in general, satisfactorily by the "exact" theoretical predic- 
tions, unlike the calculations based on the plane-wave 
impulse approximation (Ku-61) and the modified simple impulse 
approximation (Br-71, Ma-70, Mc-71) which failed to predict 
the amplitude and shape, and the amplitude of the cross 
section peak respectively. However, unlike the study of f.s.i., 
many significant discrepancies were observed between the 
"exact" theory prediction and the experiment in the p-p q.f.s. 
For example, Klein et al., who examined the D(p,2p)n break up 
at 16 keV incident proton energy and several experimental 
geometries where both the f.s.i. peaks as well as the q.f.s. 
peak appeared on the same locus curve, found a systematic 
overprediction up to 30%, of the q.f.s. peak magnitude by 
the Ebenhoh code. 

A whole summary of disagreements observed between 
the theory and experiments in the D(p,2p)n q.f.s. is given 
ine thespapersofeAnzeltongetyal .«(An-/2)%s- oi iss founds yin 
their work, that this disagreement (up to 80%) is a function 
of,incidenti energy: and. the C.M. angle..,.The, largest discre- 
pancy is found at 23 MeV bombarding energy. 

Theoretical studies of Brayshaw (Br-74) and a 
thorough analysis of Petersen et al. (Pe-74) find that q.f.s. 
shows#veryi littlessensitivity to-off-shell\ properties of the 


two-nucleon potential. Consequently, speculations were made 
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(Br-74, An-72, K1-73, Ma-71) that the existing disagreement 
is due to inadequate subtraction of the Coulomb effect in the 
p-p interaction and due to the lack of inclusion of higher 
partial waves L > 0 (only s wave potentials included) in the 
theore ti cale calcul ati.ons.. 

The unanimous conclusion of all these studies is 
the call for an accurate D(n,2n)p q.f.s. experiment where the 
Coulomb effect ambiguities are absent. Suggested are n-n 
q.-f.S. experiments at energies and geometries of the largest 
observed p-p q.f.s. disagreements (An-72). Studies of 
assymmetric q.f.s. scattering angles are suggested as 
particularly sensitive to the higher partial wave effects 
(Ma-71). 

Since the q.f.s. was found to be insensitive to off- 
shell effects it can be used to determine the on-shell 
parameters of the nucleon-nucleon potential, i.e., the 
scattering length and the effective range. An interesting 
analysis in this direction was made by Vranic et al. (Vr-/74) 
who found that, while the sensitivity of q.f.s. to the 
scattering length is poor, it is very sensitive to the 
effective range variations. Since this situation is reversed 
in the f.S.1+ experiments, a conclusion could be drawn, that 
ann value should be deduced from f.s.i. experiments while 


the effective range Prat could best be obtained from q.f.s. 


experiments. 
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Only two D(n,2n)p kinematically complete q.f.s. 
experiments have been made up to date (S1-71, Fo-/4). Both 
of these experiments were performed at 14 MeV incident 
energy and 30°, -30° q.f.s. geometry. The latter, (Fo-74) 
added a 40°, -40° scattering geometry. Unfortunately, due 
to the difficulties involved in the n-n quasi-free scattering 
experiments (see Part B), both studies contain large 
uncertainties, namely a statistical error of ~~ 30%, 


disallowing any significant conclusions. 


In view of the experience summarized in this 
chapter we decided to study neutron neutron quasi-free 
scattering in a kinematically complete D(n,2n)p experiment 
at the incident neutron energy close to the "controversial" 
23 MeV (the reason for choosing 21.5 MeV incident energy 
is explained in Chapter B.2.1). To study the effect of the 
higher partial waves, not included in the “exact" theory, 
we proposed to obtain an angular distribution of the q.f.s. 


cross section with the inclusion of assymmetric angle pairs. 


yeep ststgnos yf foo ttementd ‘ate 
“ito s(R%-07 , PC212) siteb 03 qu Sbsn 
dnabtont Vsh fi fs bomyottsq ss gx s20ns 
(Wot) ,tedsse! sit .vrtsmosp 2. +.p apie 5 8 : 
sub eyledanus otal UAT mo Se pntvetiso2 *08- oe my 
putivstiso2 sevt-fesup n-n grit ve bev foyat opis tusfttth eds of 
epyst: ‘mtesneo estbui2 dted .(4 41649 992) ednaminegxe 

S0f v to vowva faotsetist2 » \fsmen < 2otan oar! 
_enoieutanes tnooTt enpre Vos porwr sath 


aes 


otdd nt bes tvemmue sonofiaqxe ons to watv wl 
ss1t-fesup cortuen moviven Wbyse o} bebissb ional 


tnemineqxes a(nS,a)d stelqmoo yi fsotismants 6 nt ern 


Hy oe 


“Teroqavortaoa” si¢ of szald yorsas noviuen Ingsbtoat 
ypyens tysbtant Vem @.P8 pateoods not woessy ait) Vem es 
odd to da9ete and ybute of .( 1. S08 seametear aT bentstqxs ef : 


~yrosdd “So5%9" oft ni bebulont ton . 28V6W tateing wogtd 


2.%.p ond to not sudtat2th ~stupes as nieddo of | q im 
oie coe 
,evteq aslens ofisommy2es To ener ot ule i; 


‘= ide of ; 


Bs, EXPERIMENTAL S71 STEM 


B.1 INTRODUCTION 

The experimental setup used in the kinematically 
complete study of quasi-free scattering is shown in 
Fig. B.1.1. 21.5 MeV neutrons were used to bombard an array 
of three deuterated liquid scintillators (C6D6) which served 
as deuterium targets. The break-up proton from the D(n,2n)p 
reaction was stopped in the scintillator while the two 
outgoing neutrons were detected in two NE213 neutron detec- 
tors positioned at various angles with respect to the 
incident neutron direction. 

Each of the three C6D6 detectors made a different 
pair of scattering angles with respect to the two shared 
neutron detectors. Thus three independent experimental 
geometries were run simultaneously in each run (utilizing 
more effectively the bombarding neutron flux). Runs were 
taken including assymetric QFS geometries (as shown in 
Fig. B.1.1) and later sets of symmetric geometries (see 
UK fe ane eat eae a 

A description of the individual components of the 
experimental setup as well as a more detailed description 
of the experimental geometry will be given in the next two 


chapters. 
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Figure B.1.1 Schematic diagram of the experimental set-up 
in Run #1 with indication of the eight 
parameters processed by the electronics in 
eight-fold coincidence. 
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Figure B.1.2 A photograph of the symmetric quasi-free 


scattering geometry set up in Run #2. 


Figure B.3.1.2 A photograph of the three C6D6 detectors 


(right center) and the two neutron 


detectors (left) mounted in their stands. 
(Run #2). 
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Bee ne erAoT NEUTRON PRODUCTION FACILI LY 

A liquid nitrogen cooled tritium gas target 
facvlity, capable of producing 47x Tae sec fluxes of 
23 MeV neutrons from the 3H(d,n) “He reaction has been 
designed and comissioned as part of this project. The 
design, which will be described here was made to optimize 
the performance, safety, economy and reliability under the 
conditions existing in the Nuclear Physics Laboratory of 
the University of Alberta. 

The demand for high energy (23 MeV) polarized 
neutrons, raised by proposals for a kinematically complete 
study of fast neutron induced break up of the Deuteron, 
as well as polarized fast neutron scattering experiments, 
uSing the University of Alberta 7 MeV Van de Graaff could 
only be met with the (dn) He reaction. The high positive 
Q-value of this reaction (Q = +17 MeV) yields the required 


neutron energy (23 MeV with the use of a 6 MeV deuteron 


beam). This reaction, furthermore, gives high neutron 
pomerrzaclon CP O50 at 30> and P= "= 0.50 cat 90" 4esee 
Wa-70). 


beac Ves cCrIpLron-or the oi Target and Cryostat 
The relatively small thickness of the tritium 
tardet Vabopt. 1.5 mg/cm), dictated by the low energy loss 


required for the low energy deuterons (5 to 6 MeV), 
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Cross section of a cryostat designed to 
cool a 1.5 mg/cm? tritium gas target to 
liquid nitrogen temperature and to provide 
safety against tritium gas leakage by 
means of a liquid helium cryopump. 
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eliminated the possibility of using the tritium target in 


forms other than aS a gas. The physical thickness of : 


H 
solid or liquid corresponding to 1.5 mg/cm* would be only 
about 0.07 mm which would be extremely difficult to maintain. 
However, the lower the temperature of the gas the lower its 
required pressure and the smaller the required length of the 
cell. The advantages of such a situation are obvious: lower 
gas pressure means increased safety against gas cell rupture; 
shorter target brings the target geometry closer to a point 
source. Ine’ compact size of the target as particularly 
important in the polarized neutron experiments where fast 
neutrons leaving the cell at 30° to the deuteron beam 
direction are used. The shorter passage of the slow deu- 
terons through the gas also reduces the timing spread in the 
pulsed deuteron beam experiments. The 34 target of the 

FNPF consists of a cylindrical gas cell maintained at the 
liquid nitrogen temperature. The low temperature reduces the 
dimensiows sof whe cell. to i.5 cm diay x tao cm length sang 

the gas pressure to just below 1 ata to produce the required 
15 mg/cm* target thickness.. To reduce ane neat load on “tite 
cell, the deuteron beam is entered and also exited from the 
gas cell through 0.00525 mm Molybdenum windows. The deuter- 
ons were stopped on a gold lined beam stop beyond the cell 
(see Pig. Bec. l.2)) depositing only © Ol7 Woof neatiny the 


target. This arrangement allows a very efficient metal 
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Figure B.2.1.2 Detail of the tritium gas cell housing. 
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connection of the gas cell directly to a liquid nitrogen 
bath while the beam current is collected on an insulated end 
cup maintained at room temperature by an air jet on the 
outside. The body of the cell, made from 0.25 mm wall 
thickness stainless steel tube, consists of two halves 
bonded together after the two precisely punched Moly windows 
are epoxied into each half-cell from the inside and each half 
individually tested for leaks. The epoxy used in the con- 
Struction of the cell is a high conductivity silver loaded 
epoxy (EPOTEC E-32). The back half of the cell has a 1.6 mm 
0O.D. x 0.8 mm 1I.D. stainless steel capillary welded into the 
wall (see Fig. B.2.1.2) providing an uninterrupted °F gas 
filling tube between the gas cell and the tritium storage 
and handling manifold (see Chapter B.2.2.1). 

The entire gas cell is enclosed in the vacuum of 
a cryostat in which the liquid nitrogen bath, extended as 
denedtsink= to thes gas cell, Tunctions as a radiation snveld 
around a liquid Helium cryopump (see Fig. B.2.1.1). The 
presence of the cryopump provides a high degree of safety 
against tritium gas spreading beyond the volume of the 
cryostat in case of its leakage from the gas target. The 
cryopumping Surfaces cohsist of series of ‘conical’ copper 
baffles soldered to the bottom of the liquid Helium reservoir 
and are strategically positioned between the gas cell and 


the beam entrance to the cryostat. Any tritium, escaping 
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from the gas cell, would have to go around these surfaces 
before entering the beam line. The combined action of the 
baffles, the cryopumping and the large expansion volume of 
the cryostat (502%) are calculated to shatter, stop and 
freeze the shock wave of Tritium resulting from possible 

gas cell window failure. Series of tantalum collimators and 
gold masks screen the cryopumping surfaces and all constric- 
tions in the beam passage, between the entrance of the 
cryostat and the target, to avoid unwanted heat loads and 
background. The collimators are converging toward the 
target in keeping with the 3° horizontal focussing angle of 
the Mobley pulse compression system. The epoxied edges of 
the gas cell windows are masked with gold collimators, one 
being placed inside the cell (see Fig. B.2.1.2). 

The cryostat is mounted on an independent structure 
which also supported the local shielding of water, lead and 
parafin. Mallory 2000 heavy metal collimators were placed 
around the neutron exit ports. The structure provides 
alignment of the cryostat on the beam line and contains a 
self explaining control panel for the sealed vacuum system, 
tritium manifold controls and automatic liquid nitrogen 
refilling controls. The beamline immediately in front of 
the cryostat contains a 400%/sec ion pump and a system of 
vacuum baffles providing efficient differential pumping 


between the relatively poor vacuum ee Torr in the beam 
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line and the high vacuum in the cryostat » (lie 


Torr (see 
diagram in Appendix 1). This reduces both, the possibility 
of spreading the tritium into the beamline, as well as 
collection of mercury and hydrocarbons from the beamline 

on the cryopumping surfaces. The safety system, monitored 
by the interlock panel, and the fast acting valve are 
described in Appendix 1. 

The cryogenic gas target facility had undergone a 
three year period of development and tests with hydrogen 
and deuterium gas before the tritium gas was finally 
introduced into the system. The tests included simulated 
gas cell ruptures and complete gas cell filling cycles. The 
functions of the cryopump were found excellent. The liquid 
He consumption rate with 1.0 HA of beam, incident on the 
gas cell, is about 0.252%/hr, in good agreement with the 
design calculations. The current gas cell has been found 
quite rugged. It has been bombarded by 1 HA of 5 MeV 
deuteron beam for approximately 2160 hours to date, running 
for weeks at a time. The target was subjected to beam 
currents Of Up to 4 wA and pressures up to 4 ata for shore 
periods of time. Preliminary runs, to look at the neutron 
energy spectra as a function of incident deuteron energy, 
showed, that up to 5 MeV, most neutrons observed were due to 
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d + “H reactions. However, above Ed = 5 MeV the reactions 
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U0 he sent 3 


H 
a+ 34, >n +n +t 3116 


Start to produce copious amounts of low energy neutrons, 
in addition to those from d + 34 ilies THe, THis 1S. not 
a problem in elastic scattering experiments as detector 
biases are usually set well above these low energy neutrons. 
In our investigation of nN + d->Mn + n + p quasi-free 
Scattering, however, it was necessary to have a very low 
Dias tn the C6D6 "target" scintillators. Therefore, to 
lower the countrate in these detectors, the deuteron beam 
energy had to be kept below 5 MeV. With 1 ata of 31 gas in 
thesce l= at 7/7 °K the target thickness 1s anout 1.53 mg/cm<. 
This corresponds to about 80 keV energy loss for the 5 MeV 
deuterons in the gas up to the middle of the cell. With the 
inclusion of 300 keV beam energy reduction in the first Moly 
foil of the cell the reaction produced 21.5 MeV neutrons at 
zero degrees. A spectrum of 21.5 MeV neutrons obtained at 
0° with the 5 MeV deuteron beam is shown in Fig. B.2.1.3. 
The yield at zero degrees from the 3H(d,n) “He 
reaction is 4 x 10 neutrons/sr.sec. At a laboratory angle 
of 30 degrees the yield decreases to about 60% of this value, 


and with Ed = 6.5 MeV a polarization of the neutron beam 


a = 0.50 is obtained. 
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Figure B.2.1.3 Time-of-flight spectrum of direct neutrons 
at 0° produced by 6 MeV pulsed deuteron 
beam. 


mats rs Peet ard he 


‘f 
i - a a on 2 


_ 


Loan ty | | i 
Vind a 7 a 


r 


7h se mee 


41 


A set of ehergy slits, built into the cryostat 
JuUSt=an tront OF the? Gas cell (Fig. B.2.lec), With a teed= 
back to the accelerator voltage control, proved very efficient 


in holding the beam steady on the target. 


B.2.2 Tritium Storage and Filling Manifold 

A tritium storage and filling manifold was 
designed for the purpose of safe filling of the target with 
tritium gas to the needed pressure for the experimental run, 
and safe storage of the gas between experiments. 

The tritium is normally contained in a storage 
furnace, in deoxydized uranium powder, to which the tritium 
is bound chemically, forming uranium tritide UT. THES 
compound, stable at room temperature, breaks up completely 
at temperatures around 450°C (Ka-51, Ma-71). Thus heating 
the powder in the furnace to 450°C evolves the available 
tritium gas, on cooling the powder, the tritium is re- 
absorbed. This simple principle allows a pumpless closed 
connection between the storage furnace and the gas target. 

The tritium manifold consists of two such 
furnaces (Fig. B.2.2.1). These are the working furnace #1 
and the backup furnace #2 (for added absorption speed in 
case of failure) with individual seal-off valves #1 and #2 
respectively. As the activated uranium can react very 
violently with oxygen and even nitrogen, and also to enable 


purification of the tritium, the two storage furnaces are 
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Separated from the rest of the manifold by a palladium 
leak which, when heated, lets only hydrogen isotopes diffuse 
through. There is also a provision for the bypass of the 
leak by opening valve #3 and the bypass valve #4 (Fig. 
BE 22 lh 

The manifold contains three additional branches 
with their valves: 
Valve #5 - gas target port with a Bourdon combined range 

vacuum-pressure gauge. 
Valve #6 - connection to a vacuum pump and a leak detector. 
Valve #3 - fresh gas supply port to load the manifold and, 

in connection with valve #4, a bypass of the 


palladium leak. 


There are two heater elements; one for the uranium furnace 
and a second for the palladium leak. The furnace heater can 
be moved from one furnace to the other, as needed. Both 
furnaces and the palladium leak have thermocouples to deter- 
mine the temperature. 

Uranium powder is held in the ends of the inner 
1/4" bores of the stainless steel furnaces by pyrex wool 
preventing the particles from travelling beyond the reach of 
the furnace heater. The wool also efficiently fills the 
dead space of the furnace, thus reducing the amount of 
tritium, necessary in the system to fill the gas cell to the 


required pressure. 
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The entire inner volume of the manifold was 
minimized to limit the necessary amount of tritium. Most 
inner passages have 1.6 mm inner diameter while the filling 
capillary, connecting the manifold to the gas cell, has 
0.8 mm inner bore. This does not seriously impede the gas 
transfer processes as compared to the relatively slow rate 
of absorption and release of the gas. All vessels, which 
contain tritium in any stage of tritium transfer between 
the storage furnaces and the gas cell, including the 
connection from the manifold to the gas cell (plastic tube 
sheath) are double walled (and ports double closed) with the 
intermediate space either evacuated or vented and monitored 
for tritium leaks. The entire manifold is enclosed in its 
Own compact fume hood (its outer stainless steel case). 
This is evacuated and tritium leak monitored. The hand 
knobs of the six all stainless steel welded bellows valves 
are protruding outside the case (vacuum sealed through the 
bottom) for ease and simplicity of operation (see Fig. 
Bie cree 

The manifold docks into a holding mechanism at the 
bottom of the cryostat (see Fig. B.2.2.3) and can be moved 
with the cryostat as one unit off the beamline, without 
breaking the filling tube connection, or can be removed 


and stored separately for other operations (So-73). 
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Figure B.2.2.2 


Side view of the FNPF cryostat off beam- 
line with the tritium storage and handling 
manifold mounted at the bottom of the 
cryostat and connected to the gas cell 
inside the cryostat by a continuous stain- 
pea cna capillary (sheathed with plastic 
tube). 


Bottom view of the tritium storage and 
handling manifold. 
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The tritium furnaces had been flushed with 
hydrogen gas (many times repeated absorption and release 
of hydrogen from the uranium powder) before tritium was 
loaded into furnace #1. This process is necessary to clean 
the surface of active uranium and improve the gas release 
and absorption rates. 

The use of the palladium leak had to be terminated 
after many cycles due to its gradual contamination and 
eventual impermeability. This effect is well known and 
described in literature (Da-64). Allegedly, it could be 
avoided by the use of 23% silver-palladium alloy instead of 
pure palladium. The palladium leak bypass (valves #3 and 
#4) is presently being used in the tritium transfer 
processes. 

The system has been used for at least 100 fillings 
of the target to date and found very satisfactory. The 
time required to fill the cell to 1 ata is about one hour. 


The reabsorption time is only ~ 15 min. 
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B.3 DETECTORS AND GEOMETRY 
Two kinds of detectors were used in the kinemati- 
cally complete study of the D(n,2n)p breakup. 
(1) C6D6 liquid scintillators (NE230) were used as deuterium 
targets. 
(2) NE213 liquid scintillators were used to detect neutrons 


following the breakup. 


Bi Swile: Sttudye of? thes CoD6eScintidlators’ Pudser Shape 
Discrimination 


As mentioned earlier, five kinematical variables 
fully define a kinematically complete three-body breakup 
event. Thus, knowing the incident neutron energy, the 
directions of the two breakup neutrons and the energy of one 
of them should be sufficient for a kinematically complete 
D(n,2n)p experiment. However, the count rate in this type 
of experiment, being proportional to the product of two 
nuclear crosssections (7HCd.n) He and D(n,2n)p) is very low 
and random coincidences are a serious problem. Overspeci- 
ficiation of the three-body event by collecting additional 
parameters then becomes useful in separation of true events 
from the random background. 

Use of the deuterated scinillator as the deuterium 
target allowed detecting of the recoil proton and obtaining 


its recoil energy and a fast timing signal for the neutron 


time-of-flight measurement. 
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The two deuterated liquid scintillators available 
today are NE230 on deuterated benzene basis (C6D6) and NE232 
on deuterated cyclohexane basis (C6D12). Earlier studies 
reported excellent pulse shape discrimination properties of 
NE230 (Bo-72) and a higher light output than NE232 scinilla- 
tor. Although, the latter has a higher ratio of deuterons 
to carbon atoms in a unit volume, the importance of higher 
light output and better pulse shape discrimination resulted 
in the choice of NE230 scintillators as targets. 

A thorough study of the pulse shape properties of 
the C6D6 scintillator was made (So-74) prior to the proposed 
deuteron breakup experiments to investigate the possibility 
of discriminating different neutron induced interactions in 
the liquid in addition to the basic neutron-gamma ray 
discrimination. 

Liquid organic scintillators are normally used to 
allow the identifications of neutrons and gamma rays. A 
charged particle, produced in the scintillator by the neutron 
Or gamma ray causes emission of light quanta by the atoms 
along its path by means of excitation and ionization of their 
electron levels. The light produced is converted into a 
stream of electrons on a photocathode, optically coupled to 
the scintillator, and multiplied into an electronic pulse by 
a photomultiplier. The amplitude of the produced light pulse, 


in relation to the energy of the charged particle, and its 
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time characteristics depends on the history of its production 
in the scintillator. A light charged particle, like elec- 
trons, mostly excites the electrons bound in the surrounding 
molecules and their prompt deexcitation produces light with 
the wavelength in the range of sensitivity of the photo- 
cathode. Heavier nuclear particles, on the other hand, 
produce a great deal of complete ionization of the surround- 
ing molecules and the light conversion on the photocathode 
depends correspondingly on the effect of wavelength shifting 
additives in the scintillator. These complex organic 
molecules break up the more energetic ultraviolet quanta, 
produced by recombination of ions, into several softer 
("visible to photocathode") quanta through their own 
decomposition and recombination. Such light response arrives 
relatively later and its production is less efficient, as 
part of the energy is always lost in the wavelength conver- 
Sion process. The proportion of the slower and less 
efficient component of the light response is thus higher, 
the heavier and slower the particle and the higher its 
density of ionization. The charged particles, available 

for thastonvérsioh of a neutron or a gamma ray interaction 
into a light pulse in most liquid organic scintillators, are 
given by the chemical composition of the scintillator. 
Consisting almost exclusively of Hydrogen and carbon, their 
interaction with the gamma ray is by the Compton effect 


(electron) and the neutron interaction is mostly by n-p 
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elastic scattering with only small contributions from n-C 
nuclear processes up to about 20 MeV neutron energy. The 
heavily ionizing Carbon is inefficient in production of 
light and the energy of its lighter breakup fragments like 
a particles is reduced by the Q-value and by the multiple 
sharing of the kinetic energy. 

The proton is replaced by a deuteron in the 
deuterated scintillator and thus two light producing parti- 
cles detect the neutron interaction. In addition to the 
recoil deuteron from n-d elastic scattering it is also 
the proton, produced in the inelastic n-d breakup. 

The decay time of the fast component of the light 
response in the liquid organic scintillators is about 
6 nsec while the slow part typically decays in 200-400 nsec. 
For recoil protons produced in the scintillator by neutrons 
the amplitude of the slow component is about 2% of the fast 
components, whereas for electrons produced by gamma rays 
the slow component is only about 1% as large as the fast 
component. The corresponding subtle difference in the shape 
of the electronic pulse derived from such superposition of 
the fast and slow components is magnified by a suitable 
pulse shape discrimination circuit and used to distinguish 
the gamma and neutron induced events in aie S¢€nntubhator. 

Several pulse shape discrimination circuits were 
tested with the C6D6 scintillator (So-74). The selection 
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dia. NE230 scintillator was mounted on an RCA8575 photo- 
multiplier tube coupled to an ORTEC270 tube base. The zero- 
crossing time of the doubly differentiated linear signal, 
with respect to the delayed fast timing output of the 
detector as the reference time, was converted into an 
analog pulse in a time-to-amplitude converter (TAC). Its 
amplitude is the measure of the rise time of the linear 
pulse and is called Pulse Shape signal in this work. At the 
Same time, the singly differentiated linear pulse, the 
amplitude of which is the measure of the charged particle 
recoil energy in the scintillator, was collected in coinci- 
dence with the pulse shape. The two parameters were dis- 
played as a 64 x 64 channel two-dimensional histogram by 
the Honeywell 516 on-line computer. 

The tests were performed with a Pu-Be neutron 
source and the monoenergetic 23 MeV neutrons from the 
3h (d,n) He reaction. The countrate effect on the pulse 
shape discrimination was studied by adding radiation sources 
close to the detector. Reliable n-y pulse shape discrimina- 
tion was obtained up to 100 KHz. A two dimensional, 64 x 64 
channel, energy versus pulse shape plot displaying only the 
pulse shape ridge corresponding to the neutron interactions 
with the atoms of the C6D6 scintillator is displayed in 
Fig. B.3.1.1. The energy of the recoil particle runs from 


the bias of 130 keV Compton electrons at the top of the plot 
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to the maximum neutron energy of 23 MeV from the ~H(d,n 
reaction at the bottom of the picture. By moving the 
gamma ridge off scale on the left it was also possible to 
identify the splitting of the neutron interaction ridge into 
the D(n,2n)p breakup ridge (PROTON) and the deuteron ridge 
(DEUTERON), terminated by the characteristic back scattering 
peak. Further on the right is the alpha-particle ridge 
(ALPHA) from the inelastic n-C interactions. 

The experience gained in the study of the C6D6 
pulse shape properties gave confidence in the use of this 
detector to obtain two additional parameters useful in the 
examination of the n-d breakup. These are the breakup proton 
energy and the pulse shape. The proton energy and its 
direction are defined exactly in the kinematically complete 
experiment and a window can thus be set on the C6D6 energy 
parameter around the calculated value. 

Although, the pulse shape was limited to n-y 
discrimination in the quasi-free scattering region, as the 
resolution becomes poor at extremely low energies, it would 
be useful in eliminating extraneous events at other n-d 
breakup geometries where the proton energy is higher. 

For really sensitive pulse shape discrimination, 
such as needed to separate recoil protons and deuterons, 

One must use two dimensional display of energy versus pulse 


shape. Adjusting the walk on the zero-crossing discriminator 
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(TSCA) for pulse shape allows one to adjust the particle 
ridges parallel with the energy axis. Some residual 
curvature in the ridges always remains, however, especially 
at the low energy. This curvature smears up the projected 
one dimensional pulse shape spectrum and thus reduces the 
potential resolution. 

Collecting the data event by event on tape, and 
defining two dimensional curved pulse shape boundaries, as 
done in this experiment, preserves maximum sensitivity of 
the pulse shape discrimination. A periodic recording of 
two-dimensional spectra for updating the pulse shape windows 
is necessary as a protection against shifts in electronics. 
To reduce this possibility, electronic delays were replaced 
with cable delays wherever possible and periodic checks 
were made of all parameters susceptible to shifts. Cable 
delays also reduced the dead time at high counting rates. 
The use of three small detectors (3.8 cm dia. x 5 cm height) 
was also dictated by the high counting rate problem. 

The detectors were clamped in a universal stand, 
during the experiment, allowing their repositioning in all 


three space dimensions (Fig. B.3.1.2). 


B.3.2 Structure of the Double Header Neutron Detectors 
NE2V3 Tiqduird scimtrlivators were used ™ds neutron 
detectors because of their good pulse shape discrimination 


properties. Because of the low counting rate problem the 
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detectors had to be made large. The width and the height 
of the detector were chosen so that they contributed 
equally to the finite geometry energy spread in the data. 
The detectors were rectangular, 7.6 cm wide, 16.2 cm high 
and 6.35 cm thick. To improve the timing resolution of the 
detector two RCA 4522 photomultiplier tubes were mounted, 
One on the top, one on the bottom of each detector, and the 
two fast outputs were added, after passing through matching 
delavySseiNear SOQ adder (FigesBal.l)ss The: linearnsignals 
from the two ends of each detector were also added at the 
preamplifier input. 

The neutron detectors were mounted on a mobile 
stand allowing changes in angles and neutron flight distances 


for each detector individually (see Fig. B.3.1.2). 


B.3.3 The Geometry of the Experiment 

The angles, at which the neutron detectors were 
positioned in the experiment, were selected to depict the 
enhancement of the three body D(n,2n)p break up cross section 
due to the Quasi-free scattering mechanism. 

Although in the exact quasi-free scattering pro- 
CessSauinawiicn tne proton is left with zero kinetic energy. 
occurs at 41.5°, -41.5° neutron angles for an incident 
energy of 21.5 MeV, the break up is dominated by the effect 
over a wide range of angles around this geometry. Here the 


proton 1S left with non-zero share of the total kinetic 


energy. 
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In this work, the angles were limited by the need 
to obtain a detectable proton recoil in the C6D6 target 
detector. Quasi Free Scattering data from nine different 
pairs of neutron angles were collected in three separate 
runs. Three angle pairs were collected simultaneously in 
each run. 

The geometry of the first run is displayed in 
Pilgse. Belsl and B.3.3-1. ‘One symmetric (35°, =35—) “and, two 
aScvMMetrioe 43.1, =20.3°3 52.0, =21.2°j)eangle: pairs 
were studied. The neutron detectors were set at unequal 
distances from the D target. Detector #1 was set at 75 cm 
to increase the eight fold coincidence countrate, while 
detector #2 was set at 150 cm to preserve the neutron time- 
of-flight resolution with the intention of later projecting 
the data onto the neutron #2 energy axis. 

In the next two runs a symmetric geometry was used. 
The neutron detectors were set at equal distances (130 cm 
average) from the C6D6 detector, symmetrically about the 
deuteron beam axis (Fig. B.3,3.2) with.the plan of projecting 


* onto the kinematical locus. The angle pairs 


the data 
Stuavede were: raco° y =23<8- 462540. 925.00 26.0% 5 2omb? Sin 
MiemsecOndorhnedaid 26.) 6 —2o0.7 4230.07, —SORG 3 Sif. 
-31.4° in the third run. 

Attempts to keep the conditions as close as 


Dossible to the ideal Quasi Free Scattering in the first run 
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FAG Unew Dine: O.1 eA plane view of the three asymmetric 
scattering geometries running simultaneously - 
Run #1. 
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Figure B.3.3.2 A plane view of the three symmetric 
scattering geometries simultaneously 
running in Run #3. 
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resulted in many challenging experimental problems associated 
with running the C6D6 detectors at an extremely low energy 


bias. This was set at 1/4 of dl 


Am x-ray peak energy (see 
Ghapter 8.4). 

As the relative light response of the proton in 
the C6D6 scintillator to that of the electron is virtually 
unknown and unpublished for these extremely low energies it 
was discovered only after the run, when we made our own 
deuteron and proton recoil calibration measurements, that 


the 1/4 24! 


Am electron energy corresponds to ~ 170 keV 

proton energy (Fig. /B.3.3.3). At the same time, the 

minimum breakup proton energy in the three geometries 

Studied in the first run is 173 keV, 114 keV and 40 keV 

respectively. Obviously, the number of counts in the QFS 

peak area, where the proton energy drops to the minimum, is 

reduced by the effect of the bias in these three geometry 

settings. An attempt was made to correct the cross-section 

for this effect by multiplying the data with the ratio of 

the total area of the measured resolution peak of the proton 

recoil energy in the C6D6 to the part of the area above the 

detector bias. The position and the width of the resolution 

peaks were obtained by fitting a known energy deuteron recoil 

peak with the Gaussian shape. The recoil peaks were obtained 

in n-dvelastic scattering runs at various Scattering angles. 
The C6D6 energy recoil pulse was collected in a 


two dimensional histogram with the scattered neutron time- 
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Figure B.3.3.3 Light response of NE230 (C6D6) detector for 
recoiling protons and deuterons relative to 
electrons. (Compiled from: Sm-68, Pa-75, 
and Br-72 and verified by our experiment.) 
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of-flight (Fig. B.3.3.4). The two dimensional spectrum was 
truncated to include only the time-of-flight region corres- 
ponding to the n,d elastic scattering peak, background fitted, 
subtracted from the data and this projected onto the recoil 
energy axis. 

While the application of this correction can only 
be regarded as reasonable in the 35°, -35° geometry set 
all the three sets of data from the first run were valuable 
in learning about the effect of the C6D6 bias at extremely 
low proton recoils in the C6D6 scintillator. This experience 
was used in the setting up for the next runs. There the 
kinematically defined minimum proton energy was 680 keV and 
370 keV respectively. In both cases the entire finite 
resolution peak around this minimum proton energy was above 


the C6D6 detector bias setting. 
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Figure B.3.3.4 A typical shape of the deuteron recoil peak, 
in a two-dimensional plot of time-of-flight 
versus C6D6 detector pulse height, produced 
by35> elastic scattering .of?2) .50Mey 
neutrons on deuterons in the C6D6 
scintillators 


.¥5aq {fo2971 soistush sdt to saerfe Testqyt A &.€.8.8 o9uph? 


tiprt(t-Yo-smis oO Folq [snoreramrb-ows 5 nf 
beavborg .sdpret sefuq yotosteb a0d) 2uwersy 
Vom G@.f2 to pnivettiade xisesto °ae yd 

S00) sit ni enovstish no 2noituen 
voselfitafoe 


ye 


67 

B.4 ELECTRONICS 

As discussed in Chapter B.3.1, eight parameters 
were collected for each event in the kinematically complete 
D(n,2n)p experiment. The parameters in excess of the five 
variables, required for the kinematically complete defini- 
tion of the event provided additional constraints which 
were useful in the separation of the true events from the 
background. The type and origin of these parameters is 
indicated in Fig. B.1.1. They were: 
Neutron #1 time of flight - TOF 1 
Neutron #2 time of flight - TOF 2 
VEUCCOneaeCeCT0r Fl pulse wen ght 


Neutron detector #1 pulse shape 


) 

) 

) 

) 

) C6D6 detector pulse height 

) C6D6 detector pulse shape 

) Neutron detector #2 pulse height 

) Neutron detector #2 pulse shape 

All these eight parameters were processed by the electronics 
(see diagram B.4.1) in eight-fold coincidence and converted 
into eight digital words by the eight analog-to-digital 
Converters, P.0,°R, 5S. TseU and V respectively, of the 
Tennelec interface to a Honeywell 516 computer. The data 
were recorded on magnetic tape event by event and display 


monitored, at the same time, by the on-line analysis 


routine (Chapter B.5). 
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All three C6D6 detectors were treated as one by 
the electronics fanning together their fast and summing 
their slow signals. The identity of the C6D6 detector 
involved in the event was coded into the C6D6 pulse height 
parameter by means of a four-way routing (Fig. B.4.3). 

The two neutron time-of-flight spectra were 
collected as singles at the same time in an additional 
pair of analog-to-digital converters (ADC A and ADC B). 
These spectra, containing the n,d elastic scattering peak, 
were used as monitors to normalize the breakup cross-section 
using published values of the ntd differential cross-section 
(Se-72), thus eliminating the need for the exact knowledge 
of target thickness, beam current, incident neutron flux, 
CU. 

The characteristic parts of the eight fold 
coincidence circuit (Fig. B.4.1) are further described 
individually with the diagrams of the corresponding parts 
of the entire circuit diagram B.4.1. The fast timing 
signals and the slow logic signals in all of these diagrams 
are displayed as thin lines while the slow linear signals 


are indicated by thick lines. 


B.4.1 Neutron Time of Flight 
The fast anode outputs from the two headers on 
each neutron detector were time matched by nanosecond 


delays and added in a simple 502 adder. These signals then 
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were fed to an ORTEC 453 Constant Fraction Timing Discrimina- 
tor (CFTD) with the fraction set to 0.5 and the output used 
to start ORTEC 437A Time-to-Amplitude Converters (TAC). 
Delayed discriminator outputs from the ORTEC 270 tube bases 
of the C6D6 detectors were used as the stop pulse on each TAC. 
The three C6D6 detector timing signals were fanned 
together by a fast EG&G logic OR-module and later unfolded 
in the analysis by the routing bit of the C6D6 pulse height 
word. Pile up events, defined by the occurrence of fast 
outputs from two different C6D6 detectors within 500 n$ec 
resolving time, were eliminated by an anticoincidence input 


into the eight fold master gate (Fig. B.4.5). 


B.4.2 C6D6 Detector Pulse Height (Proton Recoil Energy) and 
the Four-Way Routing 


The linear outputs of the ORTEC 270 tube bases were 
fed into ORTEC 460 Delay line amplifiers. The singly 
differentiated unipolar output signals were summed by a 
summing amplifier and converted in ADCT, which set the 
routing bit of the digital word according to the combination 
of two logic signals on its routing inputs A and B. There 
are four such combinations, which was sufficient to distin- 
guish the three C6D6 detectors. The fourth routing channel 
was used for a C6H6 (normal benzene) detector in later 


experiments. 


The routing logic signals were derived from the 
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C6D6 fast signals converted into slow logic signals by a 
CFTD. The high rate of the routing signals due to extremely 
high counting rate (~ 80 kHz in each of the ee honaeeectars) 
was reduced by a coincidence requirement of at least one 
simultaneous neutron detector output. 

Four routing signals 1, 2, 3, 4 were converted into 
a combination of two routing bit A and B signals by a four- 
way routing module built in our laboratory. 

The energy bias on the C6D6 detectors was set at 


] 


1/4 fe Am source energy in the first and the third runs and 


24] 


A | A Am energy in the second run. The bias setting in 
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the neutron detectors was 1 x Co in the first run and 


1/3 ONE Compton electron energy in the second and third runs. 


B.4.3 Pulse Shape Discrimination Circuits 

The pulse shape discrimination circuit used with 
the C6D6 detectors was of the type discussed in Chapter B.3.1. 
The doubly differentiated (bipolar) output of the same 
ORTEC 460 Delay line amplifier, as used in the C6D6 pulse 
height, produced a fast timing signal in an ORTEC 450 Timing 
Single Channel Analyzer (TSCA). This was derived from its 
Zero—crossing times This’ Sign all was) used :asrethie is tartar a 
time-to-amplitude converter and the delayed C6D6 fast output 
used as the stop of the TAC. To cope with the high counting 
rate in the C6D6 detectors only cable delays (~ 1 us) were 


used on the fast C6D6 signal. 
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Further possible dead time due to a busy TAC was 
avoided by gating the TAC with the neutron detector event 
(dashed line). The DDL amplifiers were set to their 
shortest integration time (0.004 setting). The three pulse 
shape signals were summed by a summing amplifier and 
converted on ADCU. 

The neutron detector pulse shape circuit was 
identical to that of the C6D6 detectors except for the use 
of a Gate and Delay Generator (G&DG) in place of the cable 


delay, as the counting rate in these detectors was low. 


B.4.4 The Eight-Fold Master Gate 

Diagram B.4.5 indicates the logic used to define 
the eight fold coincidence gate. In addition to the double 
coincidence of the two time of flight signals the presence 
of the pulse height and pulse shape parameters from both 
neutron detectors was demanded. These two logic signals 
were obtained by application of a fully open single channel 
analyzer window on the neutron detector pulse shape TAC 
outputs. 

It should be noticed, that only the presence of the 
C6D6 fast signal was demanded, via the requirement of the 
time-of-flight TAC outputs, while the C6D6 pulse shape and 
pulse height requirements were relaxed. This allowed us to 
record on tape even the events where the proton recoil 


energy signal in the C6D6 scintillator was below the slow 
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bias, but still above the fast bias (as the fast bias was 
always set slightly below the slow one) and thus save the 
events where the C6D6 pulse shape signal was missing. This 
also saved events which would be otherwise lost due to the 
C6D6 pulse shape TAC being busy. Such events, classified as 
"zero data" by the data acquisition program (see Chapter 
B.5) were distinguished from the regular eight parameter 
events by a negative number in place of the missing word. 

As mentioned above, events having C6D6 detector 
outputs within 0.5 uw sec were cancelled by an anticoincidence 
input overlapping the three coincidence inputs in the main 
triple coincidence. This signal was derived from the double 


coincidence of two valid outputs from any two C6D6 detectors. 


B.4.5 The Monitors 

The monitor circuit is shown in Fig. B.4.6. The 
two Time of flight signals, converted into two singles TOF 
spectra via ADCA and ADCB, were obtained in the same elec- 
tronics circuit as the eight-parameter events without the 
eight fold coincidence requirement. Two features were 
different from the original circuit. An analog window was 
set on n-y pulse shape discrimination using a single channel 
analyzer on each neutron and C6D6 #1 pulse shape TAC output. 
These logic signals were used to open gates on the ADC's. 
Seconly, the C6D6 pulse shape discrimination circuit was 


duplicated for the monitors as the dynamic range of the DDA 
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amplifier, set for the low energy proton recoils, would not 
cover the elastic deuteron recoil energy in the first run 
geometry. The two different dynamic ranges were not 
necessary in the next two runs, where the proton energy was 
higher and the deuteron recoil lower at the smaller 
scattering angles. Only one C6D6 detector was used in the 


monitor circuit. 
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Figure B.4.1 Eight-parameter electronics diagram. Fast 
and logic signals indicated by thin lines, 
linear signals indicated by heavy lines. 
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Figure B.4.2 Time-of-flight circuit diagram. 
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Figure B.4.3 Electronics circuit of C6D6 pulse height 
with four way routing. 
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Figure B.4.4 Pulse shape discrimination circuits. 
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Figure B.4.5 Eight-fold coincidence gate. 
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Figure B.4.6 The electronics circuit of the monitors. 
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B.5 THE DATA ACQUISITION PROGRAM 

The General Purpose Kick Sorter software system 
(GPKS)(ref. Da-70 ), available at the University of Alberta 
Nuclear Research Centre for general on line data acquisition 
with the Honeywell 516 computer, was extended for eight- 
parameter event-by-event recording on a magnetic tape and 
modified to adapt "user" Fortran subroutines. The compiled 
program for the on line D(n,2n)p data acquisition "OLINA" 
consisted of selected GPKS routines and Fortran -subroutines 
for data collecting initialization and for a simple on-line 
data analysis. 

The block diagram of the program is shown in Fig. 
B.5.1. The teletype GPKS control dispatched the access to 
various GPKS functions. The important ones were: KSGO, 
EVIN, EVND, TAPE and DISK. EVIN opened a new event recor- 
ding file on the tape. EVND terminated the event recorded 
file by a tape mark. TAPE and DISK allowed writing of whole 
spectra regions on magnetic tape or disk as well as reenter- 
ing data from these periferals into the memory core. KSGO 
readied the ADC's waiting for an event interrupt and started 
the display of data with the data region selected by a 
letter switch. The continuous display "DO" loop was only 
exited “by a display interrupt or by an event interrupt Trom 
the sADG's2  Theddisplay interrupts, controlled by a Jight 
pen and display switches included a number of functions, 


Such as cropping and expanding of the display, summing of 
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A block diagram of the data acquisition 
program. 
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89 
peaks, etc. The displayed image could also be rotated and 
tilted - a feature particularly convenient when displaying 
two dimensional histograms. 

Two kinds of event interrupts existed. The eight- 
fold coincidence interrupt occurred when the eight-fold 
coincidence master gate logic pulse arrived at the Tennelec 
event input. A singles interrupt occurred when a monitor 
Singles time-of-flight signal arrived at ADCA or ADCB. The 
eight-fold coincidence interrupt had priority over the sin- 
gles interrupt. 

Storing of events in a buffer for the event 
recording on tape was part of the eight fold coincidence 
interrupt routine and thus it had priority over the computer 
background routines, whicn included the on line analysis. 
Two 512 word buffers, enabled by a flip-flop switch, assured 
fast recovery of the computer to accept the next event after 
writing on tape. 59 events, together with their delimiter 
words, filled one buffer. 

The on line analysis of an event, unfinished due 
to the interrupt, was completed when the computer background 


was again available. 


B.5.1 The On-Line Analysis Subroutine (ANALYZ) 
The eight fold coincidence events were accumulated 
in 32 x 32 channel two dimensional histograms of TOF 2 


versus TOF 1 by the ANALYZ subroutine. Events were sorted 
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Figure B.5.2 -A logic tree diagram of the Fortran subrouting 


'ANALYZ' contained in the data acquisition 
program as the on-line analysis of the eight- 
fold coincidence data. 
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92 
into eight such spectra, two for each C6D6 geometry. (The 
fourth C6D6 routing channel was unused.) The two regions 
for each C6D6 detector were separating into those called 
"TRUE" which satisfied a set of eight parameter windows and 
those, which did not, in regions called "RANDOM". 

The eight parameter windows were defined by light 
pen in a separate program called JSSET and then they were 
punched on a paper tape. This tape was then read in OLINA 
ied Comm Wiel zation sMODE subroutine seb ide. be Soles 

The spectra, accumulated periodically in the 
program JSSET, for the definition of the windows and their 
updating, were also written as intermittent files on the 
event recording tape. They were later used, in a more 
sophisticated way, for chronological updating of windows 
during the playback of events in the off-line analysis 
program. 

The OLINA and JSSET programs are described in 
detail in (So-76). The C6D6 pulse shape discrimination 
windows for the on-line analysis were defined, using light 
pen, as curved boundaries, intensified in the two dimensional 
pulse height vs. pulse shape display (as suggested in 


Chapterub oan.) 
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B.6 THE DATA ANALYSIS PROGRAM 

A General Off-Line Deuteron breakup Package of 
ANalysis (GOLDPAN) was written in Fortran language as part 
of this project for the off-line D(n,2n)p data reduction on 
an SDS 940 computer. 

The package consists of twelve different callable 
programs, each new program overlaying the former program 
in the memory by means of Fortran Linking procedure included 
in the MONARCH software system of the SDS 940 computer. The 
programs, linked in the GOLDPAN package, contain all rou- 
tines encountered in the process of the D(n,2n)p data reduc- 
tion. These include: definition of the eight parameter 
windows; calibration of the neutron time-of-flight and its 
conversion into meutron energy; calculation and display of 
the kinematical locus and the kinematically allowed broaden- 
ing band around it; playback of the event recorded data from 
a nine-track magnetic tape with forming and display of two- 
dimensional histograms; projections of histograms subject to 
two-dimensional cuts; storing and recall of analyzed data 
sets; plotting of the analyzed data as one-dimensional or 
two-dimensional histograms and scatter plots by a Houston 
Omnigraphic Plotter. The main interactive user control of 
the program is accomplised by a light pen and an RCT 
display. .f teletype as used tas an 1/0 device for tape, fi le 
labels and change of variable values only. An “AUTOMATIC 


PILOT" feature is included in the package, which allows tne 
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user to switch the command input to a card reader, where 
the sequence of commands for the entire analysis can be 
accumulated in a deck of punched cards. This was found 
very useful, since the analysis is usually repeated many 
time, changing only the window set and C6D6 detector numbers. 
There are seven principal links in the package, which have 
Girece Mutual access via the “JOB JUNCTION (bids. Bo.) « 
The remaining five are mere extension segments of the 
principal links, where these became too large to be con- 
tained in the directly addressable memory. They always 
return to the link from which they were called with the 
exception of Tink Ill. This returns control. to link 12° Phe 
IWnLerective control im the principal links sand )imk 11 1s 
based on a system of display pages branching into a "tree" 
from the first page in each link. Some examples of the 
disiiidyecdges are Shown In Figs. Be6w and" S.Ga3e) ebacn 
display page COnSISts OF an alphanumeric messages at the Lop 
and an eventual data display. The alphanumeric message 
contains a heading in the first line and function branch 
labels in the following rows. A light pen hit in the 
message brings another display page according to the label 
hit. A hit in the heading line returns to the previous page. 
Metight pen nit in the data display executes Tunictionse of 
the current page. These involve intensification of channel 


points or intervals of the data display in most cases. fhe 
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Figure B.6.1] The link structure of the off-line data 
reduction program. Arrows indicate the access 
between individual links. 
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Figure B.6.2 


Examples of the computer display "pages" 
included in the data reduction program. 
First row: calibration of the singles 
time-of-flight spectrum and definition of 
TOF and EY, limits. 


Second row: definition of the neutron detec- 
tor Pues shape windows. 
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Figure B.6.3 


Examples of the computer display "pages" 
included in the data reduction program. 

First row: TOF2/TOF1] histogram of "TRUE" 
events with an intensified locus (center) and 
an intensified broadening band around locus 
(right). 

Second row: TOF2/TOF1] "RANDOM" events (left), 
projection of the "RANDOM" events onto TOF] 
axis (center) and projection of "TRUE" events 
onto TOF1 axis (right). 
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System is described in detail in (So-76a) (Nuclear Research 
Center). Only the basic operations will be described here. 
The main functions of the principal links are the following. 
LINK 1 - Window Set Handler: This allows initial opening 
of the eight parameter windows by setting the lower limits 
to channel zero and the upper limits to 1023, which was the 
maximum span in all eight parameter values. The link 
facilitates writing or reading of the defined window sets 
together with all other variables of the "common" space onto 
or from a seven-track magnetic tape (unit #3). (The window 
vector is a part of the common space of each link which 
Stays undisturbed in the memory on link changing.) 

The window tape routine automatically checks the 
current position of the tape and moves the tape in the right 
direction for reading of the requested window set. It 
automatically moves the tape to its recorded end and 
assigns a new file label (window set number) when writing 
a new window set. The routines of Link 1 also include 
changing of the individual window values from the teletype. 
The last three routines always end by printing the values 
of the current window set on a line printer, including the 
window set number, for easy checks and keeping of reference. 
Link 1 does not include any data display. 
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Shape. The user branches to different display pages, 

which contain one-dimensional NE213 pulse height spectra or 
two-dimensional 32 x 64 channel pulse height (energy) versus 
pulse shape spectra (see Fig. B.6.2, second row). All of these 
sipectra™ are’ read (within the Tink) from a nine ‘track data 
tape, where they were recorded as intermittent files in a 
chronological sequence with the event-by-event recorded files 
during the experiment. The definition of a circuit is accom- 
plished by penning labels "LOWER" or "UPPER" and a subsequent 
channel hit. Both, the NE213 and C6D6 pulse shape windows 

are pulse height dependent (curved boundaries) in the off- 
line analysis, unlike the on-line analysis, where the NE213 
pulse shape limits were defined as a single window for all 
pulse heights. The current values of the limits are 

displayed as a line of intensified points during the proce- 
dure. The whole area of the data display currently enclosed 
between the "LOWER" and the "UPPER" boundaries is intensified 
following @ lirghtepen hit of Vabel*=T0TALs ~Artypical “fOTAL: 
NE213 pulse shape window is shown in Fig. B.6.2 and as the 
crosshatehnedvarea’ in Fig. Cal.Z. 

Linkeo='is the exact vanalogy of link 2efor the Co0G detectors. 
Link 4 involves both two-body and three body kinematics 
calculations to enable calibration of the time-of-flight 
Spectra and its conversion into energy - These were also 

used to define the windows on the neutron time-of-flight 


and/or on the neutron energy. 
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The program accepts input of kinematical variables 
like incident deuteron beam energy, bombarding neutron and 
the scattered neutron angles, neutron flight distances, etc. 
as well as the nsec/channel time-of-flight calibration 
obtained with a time calibrator during the experiment. The 
program facilitates comparison of calculated values with 
the shape of the singles time-of-flight spectra from ADCP 
and Q. These were recorded as intermittent files on the 
data tape during the experiment and read from the tape 
within link 4. There were five more or less pronounced peaks 
in each of these spectra, corresponding to the y-peak and 
the elastic n-d and n-c scattering peaks, each for two 
different neutron energies, 21.5 MeV and 7.9 MeV. The 7.9 
MeV group of neutrons (see Fig. B.2.1.3) resulted from the 
D(d,n)*He reaction on deuterium impurity in tne tritium 
target and on deuterium imbedded in the collimators and 
beam stop. 

Penning a certain peak of the spectrum with its 
subsequent identification by penning the corresponding 
label ,Ggoamma,,..Del 22, MeV. Ced 22 MeVsy De ly iva 9 Men. Cesl 
739 MeN. = SCen F1Gin Be G2 3 first row) caused. i nitensaticat1 on 
of; channels corresponding to the calculated position of the 
remadning, four. peaks... (hijis) gave.a werye sensitive: check «of, 
the correct calculation of the zero time-of-flight channel 
and its correct conversion into energy. This routine is 


followed by the definition of the time-of-flight and energy 
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windows in a way similar to that of links 2 and 3. 

Finally, a calculation of the three-body kinemati- 
cal locus points and its allowed limits in terms of TOF or 
energy is made, with respect to the time-of-flight or energy 
display cropping windows and stored in the "common" vector. 
Link 5 contains routines for the processing of the monitor 
Spectra recorded on the 9 track data tape after each event 
recorded run. Routines include reading of the spectra from 
tape, light pen definition of the beginning and end channels 
of the n-d elastic peak with additional points to define 
the average background on the left and on the right of the 
peak. The channel contents within the defined peak interval 
are summed up and the background is interpolated under the 
peak area. The raw sum of the peak, the net sum (raw sum 
minus background) and the sum of background inside the peak 
area are printed out, together with the statistical errors 
of the sums. At the same time a peak centroid is calculated 
and the monitor A and monitor B spectra from the individual 
runs are being added into corresponding monitor spectra for 
the whole experiment at one geometry. The centroids of the 
n-d peaks are lined up into one channel for the adding. 
Final processing of the monitor peaks was done in these 
added spectra as the statistical error there was very small. 

Links #6 and #7 differ only in the event playback 
analy sast. T8Spectica’ in = ink! 17: tarepeto nme dihian ctenmst fot iIneucron 
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verted into neutron energy and the spectra formed in terms of 

neutron 1] energy (El) and neutron 2 energy (E2). 

Eight histograms are found during the playback 
analysis in each of these two links. These are: 

(1) 64 x 64 channel two-dimensional display of E2 versus E] 
for events which satisfied the eight parameter windows 
(SE RUES) 3 

(2) 64 x 64 channel E2 versus El display of background 
events, which failed the window tests ("RANDOM"). 

(3) Spectrum of C6D6 recoil energy (pulse height) corres- 
ponding to the "TRUE" events within the allowed 
broadening band around the calculated kinematical locus. 

(4) Spectrum of C6D6 pulse heights for all the "RANDOM" 
events and for "TRUE" events lying outside the band 
around the locus. 

(5) and (7) Projections of "TRUE" events within the band 
onto the El and E2 axes respectively. 

(6) and (8) Projections of the "RANDOM" events and of "TRUE" 
events outside the band onto El and E2 axes 
respectively. 

Items. (1). (2) and (5) to (8) would read 10F in place sof ib 

in Link 7. 

These data histograms are stored in the upper 
memory of the computer (indirectly accessible 16K memory 
core) and thus they are undisturbed by changing the links 


and even when other independent programs (not using upper 
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Figure B.6.4 A logic tree of a subroutine 'ANALYZ' 
contained in link 9 of the off-line 
data reduction program. 
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SUBROUTINE ANALYZ' 
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memory) are executed. A data histogram is selected by a 
lights pens hits ine thee listtofe histogramse ony thes first display 
pages ofe links 6ior. Tinks? .to Thais) transfers rthelselected data 
region into the display buffer from the upper memory. 

Light pen hits can also intensify the calculated 
locus over the data display in histograms #1 and #2 of 
links 6 or 7 as well as intensify or change the band 
around the locus (Fig. B.6.3., first row). The bocus: and the 
band are calculated in terms of energy in link 6 and in terms 
oT time-of-flight in Vink 7. 

The whole sets of eight histograms can be stored 
On a seven track magnetic tape. Individual histograms 
from this tape can be added (or subtracted) with arbitrary 
multiplication coefficients to the corresponding histograms 
in the memory core according to the scheme: CORE = C*CORE + 
T*TAPE. The C and T are teletype input values. 

The light pen control routines also include 
projection Of the “TRUE” E2 versus El data witnin the. locus 
band onto the locus line and grouping of the projected data 
into bins of 1 MeV or other size. All data sets displayed 
in links 6/7 can also be dumped in a matrix form on a line 
printer. The “AUTOMATIC PILOT" option is selected in these 
two links. 

Link 11 contains a “HAGAR" (An-76) plotting routine modi- 
fied for plotting of the kinematical locus and/or the 


broadening band around it. The plotter offers options of 
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plotting the two-dimensional data histograms as three- 
dimensional isometric plots or as scatter plots with optional 
characters corresponding to various count number levels (see 
Chapter C - results). It also plots one-dimensional spectra 
as step-like histograms, line plots, dotted line or character 
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C. EXPERIMENTAL RESULTS 


C.] DEFINITION OF 8-PARAMETER WINDOWS 

As mentioned earlier, different C6D6 and neutron 
detector biases were set in the three different experimental 
geometry runs. The C6D6 detector biases were set at 1/4 of 


241 


the Am x-ray energy in the first run (assymmetric geometry 


anghe paiyrs 35°6°=34 3°: Heel? =2ewe°o sag? 0°=2 102°) and 
in the third run (symmetric sanglesmeo. 6% -28.6°: 30.0°, 
-30.0°; 31.4°, -31.4°) While in the second run (symmetric 
pairs 23.8°=23.. 6 Uwe 0G AS Con, suiome wsueaoeca, the C6D6 


241 24] 


bias was set at 1/3 Am energy. The 1/4 Am energy is 


equivalent to a light amplitude produced by about a 170 keV 


241 


proton GFig..B 3193). whereas the 7 3 Am bias is 


equivalent to about 220 keV proton energy. The neutron 
detector bias was set at 1 x aD Compton edge, which is 
equivalent to 3.5 MeV neutron energy, in the first run and 
ate tli/23 Nae? annihilation gamma ray Compton edge corresponding 
to about 1 MeV neutron energy in the second and third runs. 
The lower C6D6 and NE213 pulse height window limits 
were defined (by light pen in the GOLDPAN program) at the 
electronics cut-off channels corresponding to the above 
biasesy) Fig. C.l. | Shows a typical neutron detector pulse 


height versus pulse shape spectrum at the 60 C5 bias in the 


first run.  Fige Gi lee displays such? spectrum att tree l/s Nao? 
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Figure C.1.1 A two-dimensional plot of pulse height versus 
pulse shape of an NE213 detector at the lower 
bias’ of 1 x") Co. iM RUun eS 
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A two-dimensional plot of pulse height versus 
pulse shape of an NE213 detector at the lower 
bias, ofl /2 °* Nae ine RUNS 2c wand eds se tOpy 

The crosshatched area indicates a typical 
pulse shape window. A corresponding pulse: 
height spectrum of a **Na source (bottom). 
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A two-dimensional plot of pulse height versus 
pulse shape of a C6D6 (NE230) detector at the 
lower bias’ of 1/4945 °Ame inset) rouged = 7 Am 
Source - pulse height spectra in the three 
C6D6 detectors used simultaneously (under 
inset). 
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bias setting, the cross hatched area in the energy versus 
pulse shape spectrum indicating a typical neutron detector 
pulse shape window. In addition, it shows the one dimensional 
eeu pulse height spectrum in the lower part of the picture. 
Fig. C.1.3 shows a typical C6D6 pulse height versus pulse 


24] 


shape spectrum together with the 3-detector routed Am 


cella bias 


Source pulse height calibration spectra at 1/4 x 
setting. 

The time-of-flight windows were defined to crop out 
the region of interest extending from 2 MeV neutron energy in 
the first run and 1 MeV neutron energy in the second run up 
to beyond the gamma ray peak. This corresponded to approxi- 
mately a 400 channel region in each time-of-flight spectrum 
which was grouped into 64 channels. 

In the conversion of time of flight into energy the 
neutron energy limits were defined to cover a 16 MeV range 
from 2 to 18 MeV in the first run and from 1 to 17 MeV in 
the second and third runs. Thus 4 channels corresponded to a 
1 MeV bin in the E2 versus El histograms and in their 
Ey 27 axes projections. 

The neutron detector pulse shape restrictions 
together with the three detector bias settings were found 
sufficient in the rejection of background, while no apprecia- 


ble effect has been observed due to the C6D6 detector pulse 


shape restrictions. Thus these were left wide open. 
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C.2 SOURCES OF BACKGROUND IN THE EXPERIMENT 
In addition to the "TRUE" D(n,2n)p break up events 
confined to a kinematical locus in the TOF2 versus TOF] 

plOtesindteated as> curve. #le%n figs me Gnealmande Gucec. a 

number of background processes produce random coincidence 

bands and areas in the TOF2 versus TOF] or E2 versus E1 
plots. Appreciable ones are the following: 

(a) A random coincidence of an elastically scattered neutron 
from a Deuteron or Carbon nucleus, in the C6D6 detector 
(D(n,n)D or C(n,n)C), detected in one neutron detector, 
with a generally uncorrelated gamma ray (mostly) or a 
neutron from the target room background in the other 
neutron detector. While the energy or time of flight of 
the elastically scattered neutron are confined to one 
value, given by the incident energy and the scattering 
angle of the detector, the time-of-flight corresponding 
to an uncorrelated pulse in the other neutron detector 
can have any value. The result are two crossing straight 
bands parallel to the TOF or energy axes at the energy or 
TOF value corresponding to the n-D elastic scattering. 
iney are. indicated as lines #2,.93 tn EigiC cle angelic. 2, 

(b) A similar pair of crossing bands is caused by the random 
coincidence of a gamma ray, originating in the C6D6 
detector and detected in one neutron detector, with a 


generally uncorrelated pulse in the other neutron detec- 
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Figure C.2.1 Various true and background coincidence loci 
(explained in Chapter C.2) as they appear in 
a TOF2 verus’ TOFI plot. 


Figure C.2.2 Various true and background coincidence loci 
(explained in Chapter C.2) as they appear in 
an E. versus E, plot. 
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tor. =Such*bands (#45'°#5°in® Fig. C22. 1)*ontersecte the 
TOF axes at the TOF values corresponding to the gamma 
ray peak. f 
(c) A coincidence caused by the sequence of a neutron scatter- 
ing in the C6D6 detector and its C(n,n'-y) or C(n,3a0-y) 
inelastic interactions in one of the neutron detectors 
with the gamma ray produced hitting the other neutron 
detector. The same can occur with a gamma ray originating 
in the C6D6 detector and subsequently Compton scattering 
from one neutron detector to the other. As the time-of- 
flight of the gamma ray from one neutron detector to 
the other is negligible, effectively a coincidence event 
with equal time-of-flight values from both neutron 
detectors is created. In case of a Symmetric neutron 
detector geometry these events will lie on a diagonal 
line" TOFZ=TOFl or" E2=E1,° respect’ vely, {76 1in’ Figen’ Ce 22 | 
and C.2.2). More precisely, these should be two parallel 
diagonal lines separated by the gamma ray time-of-flight 
between the two neutron detectors. 
The background coincidence probability is higher 
where the bands described intersect each other causing a 
background enhancement at such locations (see Fig. C.3.2.4). 
One of these peaks was used in the analysis as an added check 
of the icorrect energy calibration.°" The calculated crossing 


position of the two n-D elastic scattering bands C7 A ne ag: 


“ as kb oe ( thas 


ti 

aie tanie noxtuay p to sonata ene ud b 
(y-0E. aya: 10 (yy! nen)o ast bas Ao 
z1039849b noyiven and Po sno At olivia oe 


novtusn reise aid patstin ‘gout yer emmep ot isiw | 
pnfisnretyo ys BM SD 6 ad tw) 1u990 ns> Smee oAT .toso9t 9b 


pitivettpo2 nod ginod Urynsundedue bas totosteb 3039 snd nt 
-to-amti ond 2A .198fto snd OF +o1oes9b novtuen sno mo1t 
ot votosteb noidven sno mort Ys% smmep ods to tdptit 
dnove sonsbfoantoa 6 yfavttostts _sfdtertpsn 2i yhio- sy 
noyjusn dtod mot 250! Vv tapi ft -te- mts Isups ddiw 


nortan aivsommye 6 Po seso nl .bstse%a 2t 2109D0g8b — 


sf 


ia! if 
: 


fedoaath 6 no sii fliw.efnavse seodt yitamost yotaetab 
L.S3. pha nr) oh) wisi tgoaqear- o[2=S3 40 fa0T=S40T oot t nie 
folisveq ows 5d bfuode s2ant .vfeetosyg avoM .(S.S.3 bas : 
Hipeft-to-swit vor simeg add vd bedersgoe zoatl femogath =~ 
.240F993.9b novsuen owt oat aaowsed ‘f ; 
tedeth 27 nivel sash fontos wie mT 7; Ai ; 7 
5 pAieyso vonte dong Soeerosnt lalla: sd ond otal w 
(ih: S28 198 sor sé) A cue : 
joan babbe 1 26 aeons 
coneey bene iwates eat 


7 


a 


121 


C.2.1, €.2.2) was intensified as a three channel corner 

together with the calculated D(n,2n)p locus in the TOF2-TOF1 

or E2/t)l display "spectra (Fig.. B.6.3). first rows. center). 
Background sources a, b and c were somewhat 

Suppressed an sthe first run as compared sto runs #2. and: 42; 

a - because of the C6D6 dynamic range excluding the elastic 

deuteron recoil, b and c - because the higher neutron detec- 

tor bias in this run wiped out lower energy secondary y- 

events and the assymmetric geometry dissolved whatever was 

left to create the diagonal band. 

(d) The above discussed bands can also occur for other 
incident neutron energy groups if present in the bom- 
barding neutron spectrum. This includes the correspond- 
ing additional d(n,2n)p loci parallel to the main breakup 
energy. locus (#8 tn Fig. C.2.1y Ci2.2). ein our case, the 
bombarding neutron spectrum contained a small impurity 
of 7.9 MeV discrete energy and a low energy continuum 
Groups. (see: Pid. Ba2. 13). 

Hhale the background sources a, b,)c can be rejected very 

efficiently by the neutron detector pulse shape restrictions 

(these events mostly involve at least one y-ray) the addition- 

al breakup background source d can not be easily rejected 

by the eight-parameter windows as these events are pure 

d(n,2n)p events, at different incident energy. However, the 

impurity of other neutron groups was almost negligible and, 


more over, the additional bands were well separated kinema- 
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tically in the E2/E1 spectrum. They were excluded in the 
projection of the data by defining a kinematically allowed 
broadening band around the main breakup locus and projecting 
the data only from this band. The corresponding projection 
band is displayed as the cross hatched area in the isometric 
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C.3 PLAYBACK OF THE EVENT RECORDED DATA WITH ANALYSIS IN 
TERMS OF THE NEUTRON TIME-OF-FLIGHT 


Typical 64 x 64 channel TOF2 versus TOF1 two- 
dimensional data histograms are presented in two forms: 
a = aS three-dimensionedeensomatriGep rows 
b - as two-dimensional scatter plots. 
Both the histograms of events accepted by the eight parameter 
window. sets... named- TRUE" or “ACCEPTED” ;-as well as—the 
rejected events ("RANDOM" or "REJECTED") are displayed for 
One geometry representing each of the three QFS runs. 

The crosshatched band in the three-dimensional 
isometric plots indicates the broadening band around the 
kinematical locus and the dotted line indicates the locus 


curve in the scatter plots. 
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Run #2. 
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CoA POPE RUINE TANALYSIS OF THE EVENT RECORDED DATA 

64 x 64 channel E2 versus E1 data histograms of 
the "TRUE" events are presented as three-dimensional isome- 
tric plots and also as two-dimensional scatter plots. One 
channel on each energy axis corresponds to 1/4 MeV interval, 
the scale covering the ranges from 2 MeV to 18 MeV in Run #1 
and fromfiZhev, to -l7) Meve im eRuns: f2eand432 

The crosshatched area in the isometric plots 
indicates the broadening band around the kinematical locus 
from which the data were projected onto either of the energy 
axes or onto a 45° straight Tine tangent to the locus at the 
OFS peak position. Thevsotirdecunve tinetie isometric plots 
and a dotted line in the scatter plot indicate the kinematical 
locus of the mean scattering geometry. Only one representa- 


tive OFS geometry is presented for each of ‘the three runs. 
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C.5 PROJECTION OF THE TWO-DIMENSIONAL DATA AND CALCULATION 
OF THE ABSOLUTE CROSS SECTION 


Gro. 1..Data erojection 

The spectra of the "TRUE" events projected onto 
the neutron energy axes from the band around the kinematical 
locus were formed automatically during the playback analysis 
of data events from the magnetic tape. 

One could, in addition, project the data within the 
band onto a 45° line, tangent to the locus curve at the QFS 
peak position, after the two-dimensional data histograms 
had been formed. The program GOLDPAN further allowed binning 
of the projected spectra channels into variable width (1 MeV 
in our case) bins together with their display and their 
Died Nt Olt. 
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The differential cross section Risnipa Glee sa the 
d(n,2n)p breakup reaction is related to the intensity 


distribution I(s 42, »&5) by the target thickness N, (cm) and 


the number of projectiles Np 
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the experimental solid angles AQ), AX, 


Thus the number of experimental counts n, in the 


B 
five-variable space interval AE, AN, AR, obtained in a run of 


duration t will be 


5 
d-o 
Peek) UN Nt ff off ae is Veet ey) 
Beale 1? e2 a qerd@ dgagisihel' 32° 2° 1 
dE, dn, de, (C-2) 


where n,(E), no(E) are the neutron detector #1 and #2 
efficiencies for a neutron of energy E. We assume the 
efficiency of the C6D6 detector for a charged particle recoil 
is 1.0 (not so in run #1 - see corrections). 

Similarly, the number N, of monitor counts, obtained 
in neutron detector 1 as singles due to the n-d elastic 
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n = 


Similarly, the number of monitor counts N, obtained when 


do, (2,) | 
Fe iariay ny (E(2)) in equation (C-3) are replaced by their 
averaged values Cas pO ves the solid angle AQ, 1s 
] ] 
NJ = N, N, t 73, ls AQ, (C-5) 


substituting for N, Ae ie AQ, from equation (C-5), into (C-4) 


we can express the average cross section in a AE. = 1 MeV 


bin around the neutron energy E. by rewriting equation (C-4) 


as 
1) n n 
Cudsalnuees te a eee eae 
i Vile 2 apg eee 
1 1 2 
The energy binning of AE. = |] MeV seemed to be a reasonable 


compromise between the variation of the differential cross 
section, and finite geometry broadening on one side and the 


statistical error on the other. 


C.5.3 Monitor Peak Summing 

The individual monitor spectra files in each run 
at.one geometry setting were added together. lining up,the 
Calculated centroids of the n-d elastic scattering peak into 
the same channel. Two such summed monitor spectra corres- 


ponding to the left and the right neutron detectors were 
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obtained in each run (#1,#2,#3) as displayed in Fig. C.5.3.1. 
These were further analyzed using the GOLDPAN routine which 
produced the net number of counts in the peak. The monitor 
peaks in runs #1 and #3 were free of any n-C elastic 


scattering contribution for all practical purposes. This 


137 


contribution was eliminated by the 1/4 C. 


C6D6 detector 
bias of the duplicated monitor circuit in run #1 and by the 


oo Mn 


combination of too small a carbon recoil energy and 1/3 
Code bias 1m run #2. 

A 23% n-C elastic contribution was observed in run 
#3. The superposition of the two peaks was reconstructed 
geometrically (see inset in Fig. C.5.3.1) on a large scale 
computer plot and the percentage of the total peak area 
corresponding to the n-d elastic scattering peak alone 
determined. Only one monitor sum is necessary in each run 
but the calculations were performed for both monitors to 
check the consistency of the results. The numbers of counts 
from the two monitors in each of the symmetric geometry runs 
were equal to within a few percent and the average value was 
ised in the Calculation of the cross section. 

The monitor spectra were taken with respect to one 
C6D6 detector only in each run (C6D6 #1 in run #2 and C6D6 #2 - 
center in runs #2 and #3). The other two sets of data in each 
run were normalized to the monitor C6D6 data set by the ratio 


of the n-d elastic peak sums in the three singles TOF spectra 
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Figure C.5.3.1 The left and right monitor spectra 
corresponding to the entire Runs #1, #2, 
and #3. The inset is showing the contribu- 
tion of n-C elastic scattering into the 
n-D elastic monitor peak and the 
reconstruction of their superposition. 
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obtained simultaneously in periodic intervals throughout each 


run. 


C.o.4 Neutron Detector Efficiency 
The neutron detector efficiency curves applied in 
the calculation of the absolute cross section, for the two 


606, eliaialeGby Sine W376, 


neutron detector bias settings of 1 x 
used In this work, are shown ih Fig. C.5.4.1. 

As no experimentally measured efficiency was yet 
available for the neutron detectors used in this experiment, 
experimentally measured and fitted efficiency curves of Drosg 
(Dr-72) were used and checked by a Monte Carlo simulation 
program (Gr-67). The Drosg data were obtained for the same 
NE213 detector thickness with the other dimensions similar 
to our detectors. The calculated efficiency curve agreed to 
within ~ 5% with the Drosg data up to the neutron energies 
where n-Carbon interactions start contributing into the light 
pulse produced in the scintillator. The computer program 
used does not contain the n-C cross sections. The upper part 
of the efficiency curve was therefore extrapolated to the 
Oe oe elie efficiency curve using the Drosg data. 

A qualitative 1)lustration of the n-Carbon contri bu- 
tion difference between the two different bias efficiency 
Curves (Fig... G.5.4.1) can be seen an Pigs. Cl land Ct toz 


as the extension of the neutron pulse shape ridge to the 


left at the low energy in Fig. ©.1.2. Whe diagonal background 
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Piqure ©.5.4.1 The neutron. detector efficiency curves 
used in this work. 


13 
Danid. Ce. Gi .@tn IFigs. €.3.2 43, C.3eee8) ko discussed iiechagter 
C.1.1, can be caused solely by the n-Carbon inelastic 
interactions in the neutron detector, because n-p interac- 
tions in the scintillator can not produce the y-ray involved 
in this type of background. Thus the diagonal band is 
further evidence for the importance of n-C interactions in a 


scintillator at the lower bias. 


Con) SCA culations “and Resa ts 

The n-d elastic scattering data of Segrave et al. 
(Se-72) were used for the n-d elastic differential cross 
section. The values obtained at 20.5 MeV and 23 MeV 
incident neutron energy were averaged to be used as the 
differential cross section at 21.5 MeV in this calculation 
since no practical difference between the two values appears 
within their error ranges. 

All data presented as tables and graphs in this 
chapter were corrected for a dead time of about 10% caused 
by the anticoincidence rejection of two simultaneous events 


occurring in two different C6D6 detectors. 
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Comparison of experimental and theoretical 
differential cross sections for the reaction 
no +deoon tnt. 1n the projecul onmoncoer 

0 heey ny 


axis in Run #1. The dashed line indicates 
the proton energy (right scale). 
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FigunesG.5. 52 


Comparison of experimental and theoretical 
differential cross sections (solid line) for 
the reaction notd > n,tnotp in projection 


onto the kinematical locus length S for the 
three angle pairs of Run #2. 
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Comparison of experimental and theoretical 
differential cross sections (solid line) 
for the reaction njtd a n,*n,*p in projec- 


tion onto the kinematical locus length S 
for the three angle pairs of Run #3. 


ARCLENGTH S (MeV) 


—— | @.=0, = 28:7° 
> 

al? 2b E4215 MeV 
: 


= 


AO” 
AQ dQ,dE 


6 Ne 12 4 &, (MeV) 


S T 
5 Ye) IS ARCLENGTH S (MeV) 


8,=8,=30:0° 


Tb Ap Eat 2ts Mev i, (7 
Ly Ny 


+ 


8 \O 2 \4 E,(MeV) 


do 
AQ AQ dE, 


ee ae 6 


AS 
d2,aN dE, 


S 4 6 8 \O 12) 4 E,CMeV) 


Hy | 4 fa, 


tobe aa 


VeoM @ig=,3 


. a I VU ai 4 a sh ae 
; a a, A tikes » aie 
/ S, ; jl eee oe 
i hd, : 


a 


SNES 06° oe 6c” 


oa 09” SHEE OP ° 
oc. ce oe Lo 
ae) iS vs 08° 
pli Ohya E Gaal vp" 
Oven Lae oy ah eo 
over ve | BL Ora 
oar €6° EV =¢ Ol’ 2 
Gat 60°¢ isn vou 
om lak Cala OCF vy 
ei Oe veel Oe 
ase 06° 60° 90° 
BC e vv Ov * 
aoa OL Cia ee 
Bd" Od - oS. 80° 


AW 4s () NeW 54s 

a a! 
u $27 unos i 

2 AONY/ °g 10) 

UA 2a. OL ————— 

abeuadry sbeuasry 


AMLANOID .0°S2- *.0°S2 GNNOYY GADYHYSAY NOIL13DIS SSOUD zZ# NNY 6 319VL 


TS Pe co cor coo io 
—- A © st to tO m~ OC 


= 
wv 
= 


= 
Ld 


igs ——- 


HUOAA 2988308 ‘worT932 ons Sk HY 


6'=@,".25* ven te 
= 21'S Ma’ ae ed 


eS 
a 


7 . ae . ne 
wae P 1 
Fa tet fe a 


oe 
tae 
. Ss 29f aie sieveles to woz tanqned > steph 
aera B>2fiordeds bone [sitesmigaqns « - y 
ad re sh NOtISG95 Bas) 99% 2064999222019 op me m= | 
apt q "OB 1+ eh “0 6 Siw Vou 2) TS, Ihe m EAvay) 
ike gots 
Bs : ape by ab reanbnF¥sax5 "gjuTooda ie 


j 
aig festsevdordebus (anit 


. 13 r4 a20tovens avitalst (s9wod)- 

. te Gfigt a nv tLe | | 
7 s 3 eoishs nojDxg « vd +2 + 10? «< 3 di 
i ar by fenyn | 


i 

é 
‘peitesb | -, Sette i” elope nottauiiorg add ‘4 : 
bey Giese adyin : 


Figure C.5.5.4 


Comparison of relevant energies and both 
experimental and theoretical differential 
cross sections for the reaction notd > n,+notp 
at EQ 2l.s MeV with > ,=0°; o5=180°, 

01 =85=25°. 

(Upper) Absolute experimental data (step 

line) and theoretical curves. 

(Lower) Relative energies En P? np? 


E for f.s.i., proton energy EY and 
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nynoc.m. (dashed line, 


the production angle 86 


right scale). 
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Figure C.5.5.5 Equivalent of Fig. C.5.5.4 for 
8,=95730°. 


Key 
A2AdAOAE, 


ENERGY (MeV) 


S) i= if eae T 
5 Xe) 1S ARCLENGTH S (MeV) 
* « Qnn Tan (fm) 
Q,=8,= 30 Be PRs vas 
is a —=— (6:4 2:86 
2 E = 2\:'S MeV ) ASB = 23-7 3h 
UY ae ee aS x coorebeo Rea 64 B34 
“y \ 
AY 
| \ 
\ 
Nig 
(e) laa T aT T i iar nah eke Tt le 7 a 
tee 4 6 B io 12 Ky  &, (MeV) 
es fe iS 20 25 


t 
‘> 


ARCLENGTH S (MeV) 


170 


\80 


i 1@) 


40 


120 


\OO 


ANGLE CDEG) 


cy vita t wf 4” 
* + .F ¥ 5 
‘ 


cat ee a ae ee hte 
(Ne) 3 SI ov A a 


Si — 


— pica “ 
sin? 


rriere (DEC) 


aly ey 


C.6 DISCUSSION OF RESULTS - CONCLUSIONS 


The averaged experimental cross sections were 
compared directly with theoretical values without taking 
into account the finite experimental resolution since this 
would have required a Monte Carlo program and an enormous 
amount of additional computer time, especially for the arc 
projection. The theoretical values (solid line) were 
obtained with the computer code of Ebenhoh (Eb-72), which 
is a separable potential approximation of the exact three- 
body theory. The theoretical cross section is projected 
onto the kinematic arc length (except for 35°, -35° data 
projected onto E2 axis), while the experimental cross 
sections are projected onto a 45° line. Deviations between 
the experimental and theoretical distribu- 
tions caused by this approximation should be very small in 
the QFS peak area since the locus curve closely coincides 
with the 45° line, except for the wings of the distribution 
where the locus has a sharp curvature. This results in a 
pileup of the projected experimental counts roughly in the 
areas where the theoretical curve predicts the minima be- 
tween the QFS peak in the center and the two njP and NoP 
Final State Interaction peaks on either end of the locus. 
Unfolding this effect would not be possible without a Monte 
Carlo program in the case of finite geometry as used in 
this experiment. The important point however is, that the 


QFS peak region remains uneffected. 
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Comparison between the experimental data obtained 
in the individual scattering angle pair geometries and the 


theoretical calculation using a = - 16.4 F and ran 2: 86 F 


nn 
for the n-n scattering length and the effective range show 
general agreement in the shape and amplitude of the QFS 
Peakywtth thermexcepti anitiof sthe 135° 235° tands2SA7 a= 428.7° 
sets. 

The discrepancy ‘observed in the 35°, -35° "set 
GhirgerGns. 521) ’could thevattributed to thie electronics 
troubles in the early stages of run #1 and to the uncertain- 
ty in the correction for the C6D6 bias interference with the 
minimum proton energy in the QFS peak (see dashed line in 
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comparison with the theoretical predictions (see various 
lites ine the top parts on fig. Cro 5es sand Cro. pols 
Although the slight tendency of the 25°, -25° set to favour 
a shorter scattering length and a longer effective range 

TS NOt convincingly sipported by the 3070" 47-3070" sset. 

a general conclusion can be drawn, that the separable 
potential approximation of the exact three body theory fits 
the data of the present experiment very well. An overall 
agreement to within about 30% both in the shape and in the 
amplitude of the data with the predictions can be 
acknowledged. 

The contributions of other possible sources of 
errors to the experimental cross section are mostly smaller 
than. the statistical error (see Table 11). 

The statistical uncertainty is lower than that 
of the only two previously existing n-n quasi-free 
scattering measurements, i.e., the data of Slaus et al. 
(S1-71) and Bovet et al. (Fo-74) obtained at the incident 
neutron energy of 14 MeV. The Slaus measurements of 30°, 
-30° quasi-free scattering was repeated in the work of 
Bovet et al. with additional measurements at 41°, -41° q.f.s. 
geometry. While the cross section at 30° agreed with the 
Slaus data, the cross section at 41° was found lower than 
the. 20° cross section bya factor of 2 tov34 slorcheck this 


discrepancy a QFS experiment is being prepared for the 
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SUMMARY OF POSSIBLE ERROR CONTRIBUTIONS IN 


THE EXPERIMENTAL CROSS SECTION 


Contribution 


Statistical’ error 


Statistical error in 
monitor sums 


Uncertainty in the 
background subtrac- 
CHOON sin mon tor 
peaks 


Error in n-d elastic 
Scattering cross 
section data 


Energy spread due to 
the finite geometry 


Uncertainty in the 
neutron detector 
absolute efficiency 


Its Effect on Cross Section 


18% and higher as displayed by 


error pars: Wh Figs... O25 17, 
CeSti2 CaS SRS oe es Rb 5 
as listed in Tables No: 2, 
A 2 oe 06. 78NCLeo ands 


negligible 


less than 5% 


~ 8% as quoted in (Se-71) 


1 MeV horizontal averaging 
bars of the experimental 
cross section data points 
as indicated by the step- 
SGlid Tanesim Fags. 9.501 
to -G.5).5 


Variations of up to ~% 20% 
found possible 


and 
35 
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41.5° - 41.5° geometry at 21.5 MeV incident neutron energy. 
A pulsed deuteron beam will be used for the production of 
the 21.5 MeV neutrons and heavy water (D,0) will be used as 
the deuterium target. 

In the meantime, QFS data have been obtained (in 
the way described in this work) also at angles (18.8°, 
=}8.8> 3 2020" 5722010" 3727 22 5N S21. 25) and’ more’ data have 
been’ taken*® at- (287 2822847 os 80. 08 Lee SU- 0c 1] 248 E514). 
The latter will hopefully allow replacing of the previous 
contradictory results at.28.7°, -28.7° geometry and thus 
give, tne possibility of obtaining a chain eof cross section 
Values: fOir 50% 0 6cOw/ seb ee. Peo ewes, Gece and 
20.0°, each an average over three adjacent sets. A least 
sauane fit to such abroad span of=se207'Sitiatasnacel= error 
data could be albe to distinguish a certain combination of 
a and Bae values best fitting the data. 

As our present results do not seem to support the 
(ip? to 80%). dis¢trepaneies* reported in’ the pp q:12s- 
experiments around 23 MeV incident energy (see Chapter 
Ne2no. 2) they could confirm the speculationsomade ihatnese 
reports that such disagreements are due to the lack of 


adequate Coulomb effect inclusion in the theoretical 


calculations. 
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An-72 


Bo-69 


Bo-72 


Br-71] 


Br-72 
Br-/4 
Ca-71] 


Ca-72 
Da-64 


Da-/0 


Dr-72 


Du=/ 1 
Eb-72 
Fa-6] 
Fo-74 


Go-70 
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APPENDIX 1. THE SAFETY SYSTEM 


In addition to the internal safety features of the 
FNPF cryostat described in Chapter B2. an external safety 
system existed (see Figure AP.1) including: Fast acting 
valve (3) triggered by a Veeco ionization vacuum gauge (2), 
system of vacuum baffles (4) and a 400 £/sec ion pump (5) 
providing efficient differential pumping, sealed vacuum 
system, vented into a duct (9) with exhaust fan at the 
roof level (10) and monitored for tritium contamination 
by a tritium monitor (11). During an experimental run 
all vacuum pumps on the beam line between the cryostat and 
the fast acting valve were outlet into a sealed storage 
tank (7), which was evacuated into the duct again after 
each run unless an accident occured. 

The following features were included in the interlock 
system: 
Vacuum gauge (2) 
Fast acting valve (3) cocked 
Mobley rotary pump outlet closed 
Fast acting valve (3) armed 
Storage tank pressure < 1 atm. (8) 
Outer case of T» manifold evacuated (16) 
T>2 monitor in stack (11) 
T, monitor of target position (12) 
Liquid He present 
Liquid Nz present 


(1 
Cryostat vacuum (17) 
Duct extract on: 


8) 


NYA ODOVUOONDAOPWND — 


ad eed od 


ond 30 2anitt net “ie ‘osha ca win 

wiatez [anv9ox9 Wb .S8 baaharaes da ero 14 

potioe tes? semtbefoat (1.98 4 yt ser ry ae a 
.(S) opuseo muyosy noteesPnot osseV 6 ed borepe ta (e) la 
(3) qmuqg sof op2\m O08 5 bee (a) 29f¥ted musey to noteye) 

mUUoBY bolsse epatamugq Ist sneer th vastots | “patbtverg 

od $6 ns? Feuerxe dttw (€) toub s ofmt “betoay ematey2 An ; 

nottentmstnos mutitas ro? bevesinom ene (ory ovat Yoon | | 

nuy Fetromiyaqxs m6 entaud .( Tf) vortaon mut tind > ae 

bas gst20N7) Sit naswtad sntl meod ant nd eqmug my uasv Hes : i- 


spsioste bafsee s otnt tofsvo stow ovisv Parry! gest aad oe 
yotts nteps toub sht osnt betsvosys. raw datdw ) ama 

‘homage jnabto258 id eeotay out fee 

Joofyednt sat at baby font orem dorusast entwor foP at 


\ ait iF Mee i 


he 


ee rie 


181 


SWANK) 


The safety system 


Figure AP.] 


AY i tai 


: 7 oe Wy bet 
77) 
: Ae i) x) 
. Kn 7 iid ne 


a 


ay . ' 
Le 
it, oF 


rc i My y 


a 


‘ 
| = on 
ms ae 


